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ABSTRACT. In this paper we show that between PDE's and crystallographic groups there is an 
unforeseen relation. In fact we prove that integral bordism groups of PDE's can be considered 
extensions of crystallographic subgroups. In this respect we can consider PDE's as extended 
crystals. Then an algebraic-topological obstruction (crystal obstruction), characterizing existence 
of global smooth solutions for smooth boundary value problems, is obtained. Applications of this 
new theory to the Ricci-flow equation and Navier-Stokes equation are given that solve some well- 
known fundamental problems. These results, are also extended to singular PDE's, introducing 
(extended crystal singular PDE's). An application to singular MHD-PDE's, is given following 
some our previous results on such equations, and showing existence of (finite times stable smooth) 
global solutions crossing critical nuclear energy production zone. 1 
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1. Introduction 

New points of view were recently introduced by us in the geometric theory of 
PDE's, by adopting some algebraic topological approaches. In particular, integral 
(co)bordism groups are seen very useful to characterize global solutions. The meth- 
ods developed by us, in the category of (non)commutative PDE's, in order to find 
integral bordism groups, allowed us to obtain, as a by-product, existence theorems 
for global solutions, in a pure geometric way. Another result that is directly related 
to the knowledge of PDE's integral bordism groups, is the possibility to character- 
ize PDE's by means of some important algebras, related to the conservation laws 
of these equations (PDE's Hopf algebras). Moreover, thanks to an algebraic char- 
acterization of PDE's, one has also a natural way to recognize quantized PDE's 
as quantum PDE's, i.e., PDE's in the category of quantum (super)manifolds, in 
the sense introduced by us in some previous works. These results have opened a 
new sector in Algebraic Topology, that we can formally define the PDE's Algebraic 
Topology. (See Refs.[46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 
63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74], and related works [2, 3, 38, 76, 77, 78].) 
Aim of the present paper is, now, to show that PDE's can be considered as extended 
crystals, in the sense that their integral bordism groups, characterizing the geomet- 
rical structure of PDE's, can be considered as extended groups of crystallographic 

^Work partially supported by Italian grants MURST "PDE's Geometry and Applications". 
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subgroups. This fundamental relation gives new insights in the PDE's geometrical 
structure understanding, and opens also new possible mathematical and physical 
interpretations of the same PDE's structure. In particular, we get a new general 
workable criterion for smooth global solutions existence satisfying smooth bound- 
ary value problems. In fact, we identify the obstruction for existence of such global 
solutions, with an algcbraic-topologic object (crystal obstruction). Since it is easy 
to handle this method in all the concrete PDE's of interest, it sheds a lot of light 
on all the PDE's theory. 

Finally we extend above results also to singular PDE's, and we recognize extended 
crystal singular PDE's. For such equations we identify algebraic-topological ob- 
structions to the existence of global (smooth) solutions solving boundary value 
problems and crossing singular points too. Applications to MHD-PDE's, as intro- 
duced in some our previous papers [69, 71], and encoding anisotropic incompressible 
nuclear plasmas dynamics are given. 

The paper, after this Introduction, has three more sections, and four appendices. 
In Section 2 we consider some fundamental mathematical properties of crystallo- 
graphic groups that will be used in Section 3. There we recall some our results 
about PDE's characterization by means of integral bordism groups. Furthermore 
we relate these groups to crystallographic groups. The main results are Theorem 
3.16, Theorem 3.18 and Theorem 3.19. The first two relate formal integrability and 
complete integrability of PDE's to crystallographic groups. (It is just this theorem 
that allows us to consider PDE's as extended crystallographic structures.) The third 
main theorem identifies an obstruction, (crystal obstruction), characterizing exis- 
tence of global smooth solutions in PDE's. Applications to some under focus PDE's 
of Ricmannian geometry (e.g., Ricci-flow equation) and Mathematical Physics (e.g., 
Navier-Stokes equation) are given that solve some well-known fundamental mathe- 
matical problems. (Further applications are given in Refs.[68, 69, 70, 71].) Section 
4 is devoted to extend above results also to singular PDE's. The main result in 
this section is Theorem 4.30 that identifies conditions in order to recognize global 
(smooth) solutions of singular PDE's crossing singular points. There we character- 
ize 0-crystal singular PDE's, i.e., singular PDE's having smooth global solutions 
crossing singular points, stable at finite times. Applications of these results to sin- 
gular MHD-PDE's, encoding anisotropic incompressible nuclear plasmas dynamics 
are given in Example 4.33. Here, by using some our previous recent results on 
MHD-PDE's, we characterize global (smooth) solutions crossing critical nuclear 
zone, i.e., where solutions pass from states without nuclear energy production, to 
states where there is nuclear energy production. The stability of such solutions is 
also considered on the ground of our recent geometric theory on PDE's stability 
and stability of global solutions of PDE's [67, 68, 69, 70, 71]. 

In appendices are collected and organized some standard informations about crys- 
tallographic groups and their subgroups in order to give to the reader a general 
map for a more easy understanding of their using in the examples considered in the 
paper. 

2. CRYSTALLOGRAPHIC GROUPS 

Definition 2.1. Let (E, (Ei,g),a) be a d- dimensional Euclidean affine space, where 
a : E x E — > E is the action mapping of the n-dimensional vector space E on the set 
E of points, and g is an Euclidean metric on E . Let us denote by A(E) =ExGL(E) 
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the affine group of E, i.e., the symmetry group of the above affine structure and 
by M(E) = E x 0(E) the group of Euclidean motions of E, i.e., the symmetry 
group of the above euclidean affine structure. (The symbol xi denotes semidircct 
product, i.e., the set is the cartesian product and the multiplication is defined as 
(a,u)(b,v) = (ab,u + av).) 2 Let us denote by R(E) = E xi 50(E) the group of all 
rigid motions of E, i.e., the symmetry group of the above oriented euclidean affine 
structure, where the orientation is the canonical one induced by the metric. One 
has the following monomorphisms of inclusions: R(E) < M(E) < A(E). One has 
the natural split exact commutative diagram 

(1) *-E > A(E) =^=s GL{E) ^1 

*- E M(E) 0(E) 1 

*- E *- R{E) 50(E) ^ 1 

The quotient groups A{E)/E = GL(E), M(E)/E ^ 0(E) and R(E)/E ^ 50(E), 
are called point groups of the corresponding groups A{E), M(E) and R(E) respec- 
tively. E is called the translations group. 

Definition 2.2. A crystallographic group is a cocompact 3 , discrete subgroup of 
the isometries of some Euclidean space. 

Definition 2.3. A d-dimensional affine crystallographic group G(d) is a subgroup 
of M(E), such that its subgroup T = G(d) HE of all pure translations is a discrete 
normal subgroup of finite index. 

If the rank of T is d, i.e., T = X d , G{d) is called a space group. The point group 
G = G(d)/T of a space group is finite, and isomorphic to a subgroup of 0(E): 
G < O(E). 

Remark 2.4. Note that the point group G of a crystallographic group G(d) does 
not necessarily can be identified with a subgroup ofG(d). In other words, G{d) is 
in general an upward extension o/T, as well as an downward extension of G, i.e., 
for a crystallographic group G(d) one has the following short exact sequence: 4 

(2) ^T— ^G(<i)^^G -1 

Thus we can give the following definition. 

Definition 2.5. We call symmorphic crystallographic groups ones G(d) such that 
the exact sequence (2) splits. 

2 Rccall that given two groups A and B and an homomorphism a : B — > Aut(A), the semidirect 
product is a group denoted by A x Q B, that is the cartesian product Ax B, with product given 
by (a, b).(a,b) = (a.a(b)(a),b.b). The semidirect product reduces to the direct product, i.e., 
A x a B = A X B = A® B, when a(b) = 1, for all b S B. In the following we will omit the symbol 
a. 

3 A discrete subgroup H C G of a topological group G, is cocompact if there is a compact subset 
K C G such that HK = G. 

^Usually one denotes such an extension simply with G(d)/T, whether we are not interested to 
emphasize the notation for G = G(d)/T. 
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Theorem 2.6 (Characterization of symmorphic crystallographic groups). The fol- 
lowing propositions are equivalent. 

(i) G(d) is a d- dimensional symmorphic crystallographic group. 

(ii) G(d) has a subgroup G < G(d) mapped by tt isomorphically onto G, i.e., G(d) = 
i{T).G andi{T)r\G = {1}. 

(iii) G(d) has a subgroup G < G(d) such that every element a G G(d) is uniquely 
expressible in the form a = i(h)g, h G T, g G G. 

(iv) The short exact sequence (2) is equivalent to the extension 

(3) *-T T xi G — — s- G >-l 

Proof. The equivalence of the propositions (i)-(iv) follows from the Definition 2.5 
and standard results of algebra. (See, e.g., [10].) □ 

Theorem 2.7 (Cohomology symmorphic crystallographic groups classes). The 
symmorphic crystallographic groups G(d), having point group G and translations 
group T, are classified by T-conjugacy classes that are in 1-1 correspondence with 
the elements ofH 1 (G;T). 

Proof. Even if this result refers to standard subjects in homological algebra, let us 
enter in some details of the proof, in order to better specify how the theorem works. 
In fact these details will be useful in the following. In the case G acts trivially on T, 
so that the group G(d) = T x G, the splitting of (2) are in 1-1 correspondence with 
homomorphisms G — > T. In general case the splitting correspond to derivations 
(crossed homomorphisms) d : G — > T satisfying d(ab) = da + a.db, for all a, b G G. 
In fact, let us consider the extension (3). A section s : G —> T xi G has the form 
s(g) = (dg,g), where d is a function G — > T. One has s(g)s(g') = (dg + g.dg' , gg'). 
So s will be a homomorphism iff d is a derivation. Two splitting si, S2, are called 
T-conjugate if there is an element h G T such that s\(g) = i(h)s2i(h)~ l , for all 
g G G. Since (h, l)(k,g)(h, 1) _1 = (h + k — gh, g) in T xi G, the conjugacy relation 
becomes dig = h + dig — gh, in terms of the derivations d\, c?2 corresponding to 
si and S2 respectively. Thus d\ and d% correspond to T-conjugate splittings iff 
their difference c?2 — d± is a function (principal derivation) G — > T of the form 
g i y gh — a for some fixed h G T. Therefore T-conjugacy classes of splittings 
of a split extension of G by T correspond to the elements of the quotient group 
Der(G, T)/P(G, T), where Der(G, T) is the abclian group of derivations G -> T, 
and P(G, T) is the group of principal derivations. On the other hand, considering 
the cochain complex G*(G,T), we see that Der(G,T) is the group of 1-cocycles 
and P(G, T) is the group of 1-coboundarics. Thus we get the theorem. □ 

Theorem 2.8 (Cohomology crystallographic group classes). The cohomological 
classification of d- dimensional crystallographic groups G(d), with point group G, is 
made by means of the first cohomology group H 1 (G;R d /Z d ). One has the natural 
isomorphism: H 1 (G;M d /Z, d ) = PP 2 (G\ r L d ). Two cohomology classes define equiva- 
lent crystallographic groups iff they are transformed one another by the normalizer 
of G in GLd(Z). Two crystallographic groups G(d), G(d)' , belong to the same 
class (arithmetical class) if their point groups, respectively G, G' , are conjugate in 
GL d (R), (m GL d (Z)). 



EXTENDED CRYSTAL PDE'S 



5 



Proof. The proof follows directly by the following standard lemmas of (co)homological 
algebra. 

Lemma 2.9 ([11]). Let G be a group and A a G-module. Let £(G,A) be the set of 
equivalence classes of extensions of G by A corresponding to the fixed action of G 
on A. Then, there is a bijection £(G,A) = H 2 (G,A). 

Lemma 2.10 ([11]). For any exact sequence 

(4) M' M M" *~ 

ofG-modules and any integer n there is a natural map S : H n (G; M") — > H n+1 (G; M') 

such that the sequence 

(5) 

H a (G; M') H°(G] M) ^ H a (G; M") — ^ H l {G; Ad') ^ ■ ■ ■ 

is exact. Furthermore if P is a projective (resp. if Q is an injective) ZG -module 
then H n (G; P) = (resp. H n (G; Q) = 0) for n> 0. 5 

In fact, it is enough to take A = T = Z d , Af = Z d , M = R d and M" = R d /Z d . □ 

Definition 2.11. We call lattice subgroups of a crystallographic group G(d) the 
set of its subgroups. 6 

Definition 2.12. We call d-dimcnsional crystal any topological space * E C E, i.e., 
contained in the d-dimensional Euclidean affine space E, having the crystallographic 
group G(d) as (extrinsic) symmetry group. We call unit cell of*E, the compact 
quotient space *E/T, having the point group as (extrinsic) symmetry group. 7 

Theorem 2.13 (Hilbert's 18th problem). For any dimension d, there are, up to 
equivalence, only finitely many- dimensional crystallographic space groups. These 
are called space-group types. We denote by [G(d)] the space-group type identified 
by the space group G(d). s 

Proof. The 18th Hilbert's problem put at the beginning of 1900, has been first 
proved in 1911 by L. Biebcrbach [8]. For modern proofs see also L. S. Charlap. 
[15]. 9 □ 

Example 2.14 (3-dimensional crystal-group types). Up to isomorphisms, there 
are 17 crystallographic groups in dimension 2 and 219 in dimension 3. However, 

^ZG is the free Z-module generated by the elements of G. The multiplication in G extends 
uniquely to a Z-bilinear product ZG X ZG — > ZG so that ZG becomes a ring called the integral 
group ring of G. A G-module A, is just a (left)ZG-modulc on the abclian group A. This action 
can be identified by a group homomorphism G — > Hom group (A). 

^The lattice subgroups of a group G is a lattice under inclusion. The identity 1 S G identifies 
the minimum {1} < G and the maximum is just G. In the following we will denote also by e = 1 
the identity of a group G. The subgroup T C G(d) is called also the Bieberbach lattice of G(d). 

7 For example if *E is identified by means of an infinite graph, (see, e.g., [89]), then the unit cell 
is the corresponding fundamental finite graph. Of course d-dimensional chains can be associated 
to such graphs, so that also the corresponding unit cells can be identified with compact d-chains. 

8 Note that a space group is characterized other than by translational and point symmetry, 
also by the metric-parameters characterizing the unit cell. Thus the number of space groups is 
necessarily infinite. 

9 See also Refs.[l, 4, 6, 7, 32, 45, 80, 84, 86], and works by M. Gromov [25] and E.Ruh [84] on 
the almost flat manifolds, that are related to such crystallographic groups. 
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if the spatial groups are considered up to conjugacy with respect to orientation- 
preserving affine transformations, their number is 230. These last can be called 
affinc space-group types. These are just the F.S. Fedorov and A. Schdnflies groups 
[21, 85]. The corresponding translation- subgroup types, or Bravais lattices, are 14- 
In Appendix A, Tab. 5, are reported the 32 affine crystallographic point-group 
types, and in Tab. 7 the 230 affine crystallographic space-group types. There are 66 
symmorphic space-group types in [G(3)]. For the other 164 G(3) cannot be identified 
with the semidirect product T x G(3)/T. It is useful, also, to know the subgroups 
corresponding to the 32 crystallographic point groups, i.e., their subgroup-lattices, 
in relation to the results of the following section where integral bordism groups of 
PDE's will be related to crystallographic subgroups. For this we report in Appendix 
B a full list for such subgroups, and in Appendix C a list of amalgamated free 
products in [G(3)]. See also in Appendix A, Tab. 8 for the 2- dimensional 17 
crystal- group types, and Appendix D for the corresponding list of subgroups. 11 

Definition 2.15 (Generalized crystallographic group). Let K be a principal ideal 
domain, H a finite group, and M a K H -module which as K -module is free of finite 
rank, and on which H acts faithfully. A generalized crystallographic group is a 
group G(d) which has a normal subgroup isomorphic to M such that the following 
conditions are satisfied: 

(i) G(d)/M = H. 

(ii) Conjugation in G(d) gives the same action of H on M. 
(hi) The extension 

(6) >■ M *-G(d) *H ^0 

does not split. We define dimension of G{d) the K-free rank d of M . We define 
holonomy group of G{d) the group H . 

Example 2.16. With K = Z one has the crystallographic groups. 

Example 2.17. A complex crystallographic group is a discrete group of affine 
transformations of a complex affine space (V, V), such that the quotient X = V/G 
is compact. This is a generalized crystallographic group with K = C. 12 

In the following we give some examples of crystallographic subgroups in dimension 
d = 3. In fact, it will be useful to know such subgroups in relation to results of the 
next section. 

Example 2.18 (The group Z®Z). This is the crystallographic subgroup, of the 
crystallographic group G(3) = Z 3 xZi, with point group C\ = Zi (triclinic syngony). 
Let us emphasize thatZQZ is crystallographic since it can be identified with the 2- 
dimensional crystallographic group G(2) = Z 2 x Zi = pi, generated by translations 
parallel to the x and y-axes, with point group %\. 

^Let A, B and C groups and B < A, B < C, the amalgamated free product AjtgC is generated 
by the elements of A and C with the common elements from B identified. 

^Let us recall that the index of a subgroup H C G, denoted [G : H], is the number of left 
cosets aH = {ah : h £ H}, (rcsp. right cosets Ha = {ha : h £ H}), of H. For a finite group G 
one has the following formula: {Lagrange's formula) [G : H] = , where o(G), rcsp. o(H), is 
the order of G, resp. H. If aH = Ha, for any a £ G, then H is said to be a normal subgroup. 
Every subgroup of index 2 is normal, and its cosets are the subgroup and its complement. 

12 See also the recent work by Bernstein and Schwarzman on the complex crystallographic 
groups [7]. 
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Example 2.19 (The groups Z 2 (J)Z„, n = 2,3,4,6). These are subgroups of the 
crystallographic groups G(3) = Z 3 x 1 n , n = 2,3,4,6, (point groups Z„, mono- 
clinic, hexagonal, tetragonal, trigonal syngony respectively) . Furth ermore, Z 2 x Z^ 7 
n = 2, 3, 4, 6 ; can be considered also subgroups of the 2 -dimensional crystallographic 
groups G(2) = Z 2 x Z n = pn, n = 2,3,4,6, (point groups 1 n , oblique, trigonal, 
square, hexagonal syngony respectively), generated by translations parallel to the x 
and y-axes, and a rotation by n/n, n — 2,3,4,6, about the origin. 

Example 2.20 (The groups Z2 @ Z2). This group coincides with the amalgamated 
free product 'Li-j^ e r Li that is generated by reflection over x-axis and reflection over 
y-axis. It is a subgroup of the crystallographic group G(3) = Z 3 xi D 4 , (point group 
-D4, tetragonal syngony). (See Appendix A and Appendix C.) The group Z2 ® Z 2 is 
also a subgroup of the 2-dimensional crystallographic groups G{2) = 1? x D4 = p4m 
and G(2) = Z 2 x D4 = p4g that have both point group D4 = Z2 x Z2, and square 
syngony. (See in Appendix A, Tab. 6 and Tab. 8.) 

Example 2.21 (The groups TL^i^D-z). (Note that D2 = Z2, see in Appendix A, 
Tab. 6). This group is generated by rotation of n, rotation by tt/2 and reflection 
over the x-axis. It can be considered a subgroup of some crystallographic group 
G(3), (see Appendix C). 

Definition 2.22. The subgroups of GL^Z) that are lattice symmetry groups are 
called Bravais subgroups. Every maximal finite subgroup of GLd(1) is a Bravais 
subgroup. 

Definition 2.23. The geometrical (arithmetical) holohedry of a crystallographic 
group G(d) is the smallest Bravais subgroup G{d) containing the point group G = 
G{d)/T ofG(d). 

Definition 2.24. A crystallographic group G(d) is said in general position if there 
is no affine transformation <fi G A(E) such that (f)G{d)f)~ l = ^G(d) is a crystallo- 
graphic group whose lattice of parallel translations has lower symmetry. 

Proposition 2.25. // the crystallographic group G(d) is in general position, then 
its holohedry G(d) is the lattice symmetry group of the parallel translations T of 
G(d). 

Definition 2.26. Two crystallographic groups belong to the same syngony, (Bravais 
type), if their geometrical (arithmetical) holohedries coincide. 

Example 2.27. For the 3- dimensional case there are 73 arithmetic crystal classes. 
The space group types with the same point group symmetry and the same type of 
centering belong to the same aritmetic crystal class. An arithmetic crystal class 
is indicated by the crystal symbol of the corresponding point group followed by the 
symbol of the lattice. Furthermore G(3) has 7 syngonies and 14 Bravais types of 
crystallographic groups. (See in Appendix A, Tab. 7. The number between brackets 
() after the symbol of the point group is the number of space-group types with that 
point group.) The geometric crystal classes with the point symmetries of the lattices 
are called holohedries and are 7. The other 25 geometric crystal classes are called 
mcrhocdries. The Bravais classes (or Bravais arithmetic crystallographic classes ) 
are the arithmetic crystal classes with the point symmetry of the lattice. The Bravais 
types of lattices and the Bravais classes have the same point symmetry. 
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3. INTEGRAL BORDISM GROUPS VS. CRYSTALLOGRAPHIC 

GROUPS 

In the following we shall relate above crystallographic groups to the geometric 
structure of PDE's. More precisely in some previous works we have characterized 
the structure of global solutions of PDE's by means of integral bordism groups. Let 
us start with the bordism groups. 13 

Theorem 3.1 (Bordism groups vs. crystallographic groups). Bordism groups of 
closed compact smooth manifolds can be considered as subgroups of crystallographic 
ones. More precisely one has the following: 

(i) To each nonoriented bordism group fl n , can be canonically associated a crys- 
tallographic group G(q), (crystal-group of f2„), for a suitable integer q, (crystal 
dimension of Ct n ), such that one has the following split short exact sequence: 

(7) >- W 5=^= G(q) == fi n ^ 

So £l n is at the same time a subgroup of its crystal- group, as well as an extension 
of this last. If there are many crystal groups satisfying condition in (7), then we 
call respectively crystal dimension and crystal group of f2„ the littlest one. 

(ii) To each oriented bordism group + f2„, n ^ mod A, can be canonically associated 
a crystallographic group G{d), (crystal-group) of + 0„, for a suitable integer d, 
crystal dimension) of + tt n , such that one has the following short exact sequence: 

(8) *-+f! n *-G(q) ^G(d)/+n n ^0 

So + £l n is a subgroup of its crystal- group. 

Proof. Let us recall the structure of bordism groups. 

Lemma 3.2 (Pontrjagin-Thom-Wall [88, 91, 96]). A closed n-dimensional smooth 
manifold V, belonging to the category of smooth differ entiable manifolds, is bordant 
in this category, i.e., V = DM , for some smooth (n+1) -dimensional manifold M , iff 

the Stiefel- Whitney numbers < uii 1 ■ • ■ u>i p , /iy > are all zero, where i\ H \-i p = n 

is any partition of n and nv is the fundamental class of V . Furthermore, the 
bordism group £l n of n-dimensional smooth manifolds is a finite abelian torsion 
group of the form: 

(9) On = Z 2 0--v®Z 2 , 

where q is the number of nondyadic partitions of n. Two smooth closed n- 
dimensional manifolds belong to the same bordism class iff all their corresponding 
Stiefel- Whitney numbers are equal. Furthermore, the bordism group + il„ of closed 
n-dimensional oriented smooth manifolds is a finitely generated abelian group of the 
form: 

(10) +n nSZ 0...0z0z 2 0...0z 2 , 



13 For general informations on bordism groups, and related problems in differential topology, 
see, e.g., Refs.[31, 39, 41, 52, 79, 83, 87, 88, 91, 92, 93, 94, 96, 97]. 

partition (ii, ■ ■ ■ , i r ) of n is nondyadic if none of the ip are of the form 2 s — 1. 
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where infinite cyclic summands can occur only if n = mod 4. Two smooth closed 
oriented n-dimensional manifolds belong to the same bordism class iff all their cor- 
responding Stief el- Whitney and Pontrjagin numbers are equal. 15 

Let us recall that a group G is cyclic iff it is generated by a single element. All 
cyclic groups G are isomorphic either to Z p , p S Z + or Z. G = Z p iff there is some 
finite integer q such that g q = e, for each g € G. Here e is the unit of G. A group G 
is virtually cyclic if it has a cyclic subgroup H of finite index. All finite groups are 
virtually cyclic, since the trivial subgroup H — {e} is cyclic. Therefore, to compute 
the finite virtually cyclic subgroups of the d-dimensional crystallographic groups 
is equivalent to compute the finite subgroups. Of course can be there also infinite 
virtually cyclic subgroups. For these we can use the following lemma. 

Lemma 3.3 (Scott & Wall [87]). Given a group G, the only infinite virtually cyclic 
subgroups of G will be semidirect products H x Q Z and amalgamated free products 
A+bC for H,A,B,C < G. 

The finite subgroups of the crystallographic groups may be derived exclusively from 
their point groups. The following lemmas are useful to gain the proof. 

Lemma 3.4. (a) Let f = (0, b) be a finite- ordered element of a crystallographic 
group G(d). Furthermore, let x = (v, 1) be a translation. Then, f commutes with 
x iff b(v) = v, i.e., b fixes v. 

(b) F X Z is a subgroup of some d-dimensional crystallographic group iff F is a 
subgroup of some (d — 1)- dimensional crystallographic group. 

Proof, (a) It follows directly by computation. 

(b) Let assume that F is a subgroup of some (d — l)-dimensional crystallographic 
group G(d — 1), then F x Z is a subgroup of the d-dimensional crystallographic 
group G(d — 1) x Z. Vice versa, let F x Z be a subgroup of some d-dimensional 
crystallographic group G'(d). Since in this direct product the generator x G Z 
commutes with all / G F and x belongs to the Bicberbach lattice T' of G'(d), 
there arc at most (d — 1) independent elements ( 6 T' which do not commute with 
a given / G F. Therefore, F must be a subgroup of some (d — l)-dimcnsional 
crystallographic group. □ 

Lemma 3.5. // a crystallographic group G(d) admits a subgroup F x Q Z for some 
finite group F and some homomorphism a : Z — > Aut(F), then G(d) also admits a 
subgroup FxZ. 

Proof. Since Z is cyclic, a is completely determined by the a(l). Since F is finite, 
Aut(F) is also finite, hence one has a(l) q = 1 for some finite q. Therefore the 
elements y = x pq G Z, Vp G Z, commute with any / G F, hence F x Q Z contains 
as a subgroup FxZ, identified with the couples (/, y), were y are above defined 
elements. □ 

Lemma 3.6. Let H = F x Q Z, Z =< x >, be a subgroup of a crystallographic 
group G(d). The elements f 6 F fix the shift vectors x q G T. 



5 Pontrjagin numbers are determined by means of homonymous characteristic classes belong- 
ing to H'(BG, Z), where BG is the classifying space for G-bundles, with G = S p (n). See, e.g., 
Refs.[39, 41, 83, 88, 91, 92, 93, 96, 97]. 
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Proof. According to Lemma 3.5 we can consider the subgroup K = F x Z < H . 
Then the elements z £ Z of K correspond to translation vectors x q £ Z < H. On 
the other hand, since K is a direct product F x Z. all / £ F commute with all 
z £ K. By Lemma 3.4(a), all / £ F fix all £ H. □ 

Lemma 3.7 (Alperin & Bell [4]). Let F and H be groups, let a : H — > Aut(F) 
b a homomorphism, and (f> £ Aut{F). If (p is the inner automorphism of Aut(F) 
induced by <p, then F x^ oq H = F x Q H . 

( This means that if we identify the conjugacy classes of authomorphisms of a given 
group, we need to consider only one element of each class to evaluate the candidacy 
of all automorphisms in that class.) 

Lemma 3.8. If the presentation of the amalgamated free product contains two or 
more elements of order two that do not commute, then the amalgamated free product 
is not a subgroup of any three dimensional crystallographic group. 

Proof. In three dimension there are only three possible elements of order two, in- 
versions, 7r rotation, and reflection. All these symmetry commute. Therefore, an 
amalgamated free product with two order two elements that do not commute cannot 
exist in 3-dimensions. □ 

Let us first note that for bordism groups identified with some finite or infinite cyclic 
groups, theorem is surely true by considering the following two standard lemmas. 

Lemma 3.9. If H is a finite subgroup of a group G, every element a £ H generates 
a finite cyclic subgroup < a >= Z„ C H , where n is the order of a, anda^ 1 = a n , 
or equivalently a n = e, where e is the unit of H (and also that of G.) 

Lemma 3.10. Every element a of a group G generates a cyclic subgroup < a >C G. 
If a has infinite order, then < a >= Z. 

Thus if Clp = 7Li. : or + Q p = Z, it follows that theorem is proved. 

Now, let us consider the more general situation. We shall consider the following 

lemma. 

Lemma 3.11. The group I s x Z| can be considered a crystallographic group in the 
Euclidean space R s . 

Proof. The Z s -conjugacy classes of splittings of the split extension 

(11) ^Z s xZ| ^Z| ^0 

are in 1-1 correspondence with the elements of 

( H^Z^Z*) = Hom- L [H 1 {Z s 2 ;Z);Z s ) 

(12) <^ £* (Hom^H^; Z);Z)) S 
[ S (ff O m z (0;Z)) s ^0. 

We have used the fact that for any finite cyclic group K of order i one has 
(13) 

i = ] ( i = 

Hi{K;Z) = { Z-2 i odd \ => H l {Z 2 \ Z) = I Z 2 i odd 

i > even. , i > even. 
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Furthermore, the Kiinneth theorem for groups allows us to write the unnatural 

isomorphism: 

(14) 

H s (GixG 2 ; Z) = fl p (G 1 ;Z)® z fl,(G 2 ;Z) T r z (i/ p (G i; Z), tf,(G 2 ; Z)) 

p+q=s p+q=s—l 

for any two groups Gi and G 2 . Taking into account that Tor 1 {A, B) = for projec- 
tive Z-module A (or £?), we get that iJ 1 (Z|; Z) = 0. Thus we conclude that there is 
an unique split in (11). Furthermore, the set of equivalence classes of s-dimensional 
crystallographic groups with such a point group are in 1-1 correspondence with 

(15) { =0 se{ i,... , r} -Hom z (Z 2 ®Z 2 ;Z?) 
[ S0 ae{1 [...] r} Fom z (Z 2 ®Z a ;Z)«, 

where r = r(g). □ 

Therefore if Vt n = Z| = Z 2 x • • - a • • • x Z 2 it can be identified with a point group 
G of a crystallographic group G(q), belonging to one of these equivalence classes, 
such that G(q) < M(M. q ). So one has the exact sequence 

(16) 1* G(q) ^ Q„ 1 

that proves that fl n admits a crystallographic group as its extension. Now, since 
G(q) contains also as subgroup Z q x Z|, that contains as subgroup Z|, it follows that 
G(q) contains also as subgroup f2„. So one has also the following exact sequence: 

(17) o — n n > G(q) »- G{q)/n n ^ 

Since G(q)/£l n = Z 9 , it follows that sequence (17) is the split sequence of (16), and 
vice versa. 

Let us consider, now, the oriented case. Let us exclude the case n = mod 4. Then 
we can write + f2„ = Z r x Z|. Let us assume r > s. Then we can consider in M. r the 
crystallographic group G(r) = Z r XI Z|. This contains as subgroup Z r x Z s 2 = + il n . 
Therefore one has the following exact sequence: 

(18) ^+Q n ^G(r) ^G(r)/+n n ^ 

Let us assume, now, that r < s, then one has the following sequence of subgroups: 

(19) + 0„ =Z r xZJ< Z s x Z s 2 < Z s x Zj^ = G(s). 
So we have the following short exact sequence 

(20) o — ^+n n *~G(a) ^G{s)/+n n ^0 

Therefore we can conclude that + fl n is a subgroup of the crystallographic group 
G(d) = Z d x Z c l with d = max{r, a}. □ 

Remark 3.12. It is important emphasize that crystallographic groups and bordism 
groups, even if related by above theorem, are in general different groups. In other 
words, it is impossible identify any crystallographic group with some bordism group, 
since the first one is in general nonabelian, instead the bordism groups are abelian 
groups. 

We can extend above proof also by including bordism groups relatively to some 
manifold. 
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Theorem 3.13. Bordism groups relative to smooth manifolds can be considered as 
extensions of crystallographic subgroups. 

Proof. Let us recall that a k-cycle of M be a couple (TV,/), where N is a k- 
dimensional closed (oriented) manifold and / : N — > M is a differentiable mapping. 
A group of cycles (TV, f) of an n-dimensional manifold M is the set of formal sums 
^2ANi, fi), where (iVj, fi) are cycles of M. The quotient of this group by the cy- 
cles equivalent to zero, i.e., the boundaries, gives the bordism groups Q S (M). We 
define relative bordisms Q S (X, Y), for any pair of manifolds (X,Y), Fcl, where 
the boundaries are constrained to belong to Y. Similarly we define the oriented 
bordism groups + fl s (M) and + Q_ S (X, Y). One has O s (*) — and + H S (*) — + ^s- 
For bordisms, the theorem of invariancc of homotopy is valid. Furthermore, for any 
CW-pair (X,Y), Y C X, one has the isomorphisms: Q^(X,Y) = Q S (X/Y), s > 0. 
One has a natural group-homomorphism £l s (X) — > H s (X;l2)- This is an isomor- 
phism for 8 = 1. In general, Q S (X) ^ H S (X;Z2). In fact one has the following 
lemma. 

Lemma 3.14 (Quillen [79]). One has the canonical isomorphism: 

(21) Qp(X)= H r (X;Z 2 )®z 2 Q s . 

r+s—p 

In particular, as Qo = li and f^i = 0, we get J7 1 (X) = H\(X; 1<i). Note that for 
contractible manifolds, H S (X) = 0. for s > 0, but fl s (X) cannot be trivial for any 
s > 0. So, in general, fl s (X) H s (X;l2). 

So we get the following short exact sequence: 

(22) o — ^K P (x) — ^n p (x) — ^n p — 

where K p (X) = 

r+s=p.r>o Hr(,X> ^2) ^s- Therefore by using Theorem 3.1 
we get the proof soon. □ 

Let us, now, consider a relation between PDE's and crystallographic groups. This 
will give us also a new classification of PDE's on the ground of their integral bordism 
groups. 

Definition 3.15. We say that a PDE E k C J*(W) is an extended 0-crystal PDE, 
if its integral bordism group is zero. 

The first main theorem is the following one relating the integrability properties of 
a PDE to crystallographic groups. 

Theorem 3.16 (Crystal structure of PDE's). Let E k C J„(W) be a formally 
integrable and completely integrable PDE, such that dim E k > 2n + 1 . Then its 
integral bordism group £1^—1 * s an extension of some crystallographic subgroup. 
Furthermore if W is contractible, then E k is an extended 0-crystal PDE. 

Proof. Let us first recall some definitions and results about integral bordism groups 
of PDE's. (For details see Refs.[52, 54, 55, 56, 58, 60, 61], and the following related 
works [2, 3].) Let W be an (n+?n)-dimcnsional smooth manifold with fiber structure 
7T : W — > M over a n-dimensional smooth manifold M. Let E k C J*(W) be 
a PDE of order k, for n-dimcnsional submanifolds of W. For an "admissible" p- 
dimensional, p £ {0, • • • , n— 1}, integral manifold N C E k , we mean a p-dimcnsional 
smooth submanifold of E k , contained in an admissible integral manifold V <Z E k , ot 
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dimension n, i.e. a solution of Ek, that can be deformed into V, in such a way that 
the deformed manifold N is diffcomorphic to its projection X = nk,o(N) C W . In 
such a case the /c-prolongation X 1 -^ = N. The existence of p-dimensional admissible 
integral manifolds N C Ef. is obtained solving Cauchy problems of dimension p G 
{0, ■ • • , n — 1}, i.e., finding n-dimcnsional admissible integral manifolds (solutions) 
of a PDE Ek C J^(W), that contain iV. Existence theorems for such solutions 
can be studied in the framework of the geometric theory of PDE's. For a modern 
approach, founded on webs structures see also [2, 3]. A geometric way to study the 
structure of global solutions of PDE's, is to consider their integral bordism groups. 
Let Ni C Ek, i = 1, 2 be two (n— l)-dimensional compact closed admissible integral 
manifolds. Then, we say that they are Ek-bordant if there exists a solution V C Ek, 
such that dV = N\ 1J N2 (where 1J denotes disjoint union). We write N ~£ fc 
This is an equivalence relation and we will denote by £l n ti the set of all E^-bordism 
classes [N]E k of (n— l)-dimensional compact closed admissible integral submanifolds 
of Ek- The operation of taking disjoint union defines a structure of abelian group on 
I n or der to distinguish between integral bordism groups where the bording 
manifolds are smooth, (resp. singular, resp. weak), we shall use also the following 
symbols (resp. s , resp. fl^ f i lw ). 16 Let us first consider the integral 

bordism group f2^*i w (or £1^-1 s)i f° r wea k solutions (or for singular solutions). 
We shall use Theorem 2.15 in [60], that we report below to be more direct. (See 
also [2, 3].) 

Theorem 2.15 in [60]. Let Ek C J^(W) be a formally and also completely inte- 
grable PDE, such that dimEk > 2n + 1. Then one has the following canonical 
isomorphism: 

(23) Q PtW ^ H r (W;Z 2 )® Z2 fl s . 

r,s,r+s— 

Furthermore, if W is an affine fiber bundle ix : W — > M over a n-dimensional 
manifold M , one has the isomorphisms: 

(24) 0^ S Qp(M) = H r (M; Z 2 ) ©z 2 fl e . 

r,s,r-\-8=p 

So we can write for the weak integral bordism group 

(25) ttn-x, w = H r (W;Z 2 )®z 2 n s . 

r+s—n— 1 



°Let us recall that weak solutions, are solutions V, where the set of singular points of V, 

contains also discontinuity points, q, q' 6 V, with TTkfiil) = 0(9') =06 W, or iTkil) = ^k(l') = 
p S M. We denote such a set by £(V)g C T,(V), and, in such cases we shall talk more precisely of 
singular boundary of V, like (dV)s = dV \ £(V) However for abuse of notation we shall denote 
{dV)s, (resp. T,(V)s), simply by (8V), (resp. S(V)), also if no confusion can arise. Solutions 
with such singular points are of great importance and must be included in a geometric theory 
of PDE's too [60]. Let us also emphasize that singular solutions can be identified with integral 
n-chains in E^, and in this category can be considered also fractal solutions, i.e., solutions with 
sectional fractal or multifractal geometry. (For fractal geometry see, e.g., [20, 35, 40].) 

^The bording solutions considered for the bordism groups are singular solutions if the symbols 
3fc and gk+i arc different from zero, and for singular-weak solutions in the general case. Here we 
have denoted fl p {X) the p-bordism group of a manifold X. For informations on such structure of 
the algebraic topology see. e.g. [31, 52, 79, 83, 88, 90, 91, 92, 93, 96, 97]. 
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Then one has that fL 



is an extension of the bordism group O n -i' Hence by 



using Theorem 3.1 we get the proof. 

Let us now consider the integral bordism group Fl^-i f° r smooth solutions. This 
bordism group is related to previous one, and to singular bordism group £1^-1 si by 
means of the exact commutative diagram (26). Furthermore, the relation between 
and f2 n _i is given by means of the exact commutative diagram (27), 



"ra-li "n— l.to 



where K, 



n—l,w 



ker(a), K*l 



this we get that also £1^-1 can be considered an extension of fl n -i if fl^-i w is 
so. Therefore we can apply Theorem 3.1 also to the integral bordism group for 
smooth solutions whether it can be applied to the integral bordism group for weak 
(or singular) solutions. Therefore theorem is proved. 



l,w;n— 1 



ker(6), K E \ 



ker(c), with c = b o a. From 



(26) 



K 



E k 

,— l,wf (s,w) 



K 



Th— 1 , S 



n 



E k 



o 



n—l,w 



K 



E k 

n—l,s,w 



ni 



E k 

n—l.w 



(27) 



n-l 



n— l,w:n— 1 u 



n— 



□ 



Definition 3.17. We say that a PDE E k C J„(W) is an extended crystal PDE, if 
conditions of above theorem are verified. We define crystal group of Ek the littlest 
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crystal group such Theorem 3.16 is satisfied. The corresponding dimension will be 
called crystal dimension of Ek ■ 

In the following we relate crystal structure of PDE's to the existence of global 
smooth solutions, identifying an algcbraic-topological obstruction. 

Theorem 3.18. Let Ek C J„(W) be a formally integrable and completely integrable 
PDE. Then, in the algebra U n _i(E k ) = Map{Q^Li, R), (Hopf algebra of E k ), there 
is a subalgebra, (crystal Hopf algebra) of Ef. . On such an algebra we can represent 
the algebra R G ( d ) associated to the crystal group G(d) of Ek- (This justifies the 
name.) We call crystal conservation laws of Ek the elements of its crystal Hopf 
algebra. 

Proof. In fact the short exact sequence 

(28) K E n 1 1 Of* l fin- 1 

obtained from the commutative diagram in (27), identifies for duality the following 
sequence 



(29) 0^ Kf*!-*— :-Hn-i(£?fc)-&— 

On the other hand from the short exact sequence 

(30) fi„-i — G(d) i G(d)/n n -i *- 1 

identifying f2„_i with a crystallographic subgroup, we get for duality the following 
sequence 

(31) -= R n -i RG(d) K G(d)/a„_! < o . 

So we can identify the crystal Hopf algebra of Ek with R"™- 1 . On such an algebra 
we can represent all the Hopf algebra R G ( d ) associated to the crystal group G(d) of 
E k . 18 □ 

Theorem 3.19. Let Ek C J^{W) be a formally integrable and completely integrable 
PDE. Then, the obstruction to find global smooth solutions of Ek can be identified 
with the quotient H„_i(i? 00 )/R f2 "- 1 . 

Proof. Let us first consider the following lemma that gives some criteria to recognize 
(n — l)-dimensional admissible integral manifolds in Ek- 

Lemma 3.20 (Cauchy problem solutions criteria). 1) Let Ek C J*(W) be a for- 
mally integrable and completely integrable PDE on the fiber bundle tt : W — > M, 
dirnH^ = m + n, dim M = n. Let N C Ek be a smooth (n — 1)- dimensional in- 
tegral manifold, diffeomorphic to X = 7tk,o{N) C W, such that X C Y, where Y 



Let us remark that sequences (29) and (31) do not necessitate to be exact, but are always 
partially exact. In fact, even if it does not necessitate that d e* =0 and /, oj, = 0, (hence 
neither im (e*) = ker(c(*) and im (g*) = ker(/*)), one has that e* and g* are monomorphisms and 
d* and /* are epimorphisms. This is enough for our proof. Let us recall also that H n — l{Ek) = 
Mapifl -jii^R) is an Hopf algebra in extended sense, i.e. it contains the Hopf algebra R^™- 1 as 
a subalgebra. (See also [55].) 
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is a smooth n- dimensional submanifold of W , satisfying the condition that its to- 
order prolongation y( fe ) c Jn~i(W), contains N, (0-order admissibility condition). 
Then, there exists a (weak, or singular, or smooth) solution V C Ek, such that 
N C V. 

2) Furthermore, if the symbols gk and gk+i, of Ek and Ek+i respectively, are dif- 
ferent from zero, then V can be a singular or smooth solution. Moreover, if there 
exists a nonzero smooth vector field £ : Ek — > TEk, transversal to N , and charac- 
teristic, at least for a sub-equation of Ek, then a smooth solution V passing through 
N can be built by means of the flow 4> associated to £. For suitable initial Cauchy 
integral manifolds, solutions can be built by using infinitesimal symmetries of the 
closed ideal encoding Ek- (For details see Theorem 2.10 in [55]. ) 

3) Finally, let Ek — > Ek-i = 7Tfc,fc-i(-Efc)> be an affine subbundle of TTk,k-l ■ 
•^n(W) — ^ J-n \W)> w tth associated vector bundle Ttk—i '■ 9k E}--\, where 
T^kk-i = fc-l °Kk; with fffc : gk — > Ek is the canonical projection. Let N C Ek be 
a smooth (n — 1) -dimensional integral manifold, diffeomorphic to X = TTk,o(N) C W 
and satisfying the 0-order admissibility condition. Then, there exists a (singular, 
or smooth) solution V C Ek, such that N C V . 

An integral manifold N C Ek, as above defined and contained into a solution V C 
Ek, is called admissible. 

Proof. Since N is an (to— l)-dimensional smooth integral manifold, diffeomorphic to 
X = TTk.o(N) C W , satisfying the 0-order admissibility condition, we can consider 
N C C J*(W), where Y is just a n-dimensional smooth manifold of W, 

containing X. In general <f_ Ek, but taking into account that Ek is formally 
integrablc and completely intcgrablc, we get that (Ek)+ r is a strong retract of 
Jt r (W), Vr > 0. Then, we can deform C Jn +r {W) into E k+r , obtaining 

a (weak, or singular, or smooth) solution Y C Ek+ r , passing for = X^ k+r \ 
Then nk +r .k(Y) = V C Ek is a solution of Ek, passing through N. In particular, 
if gk 7^ 0, and gk+r ^ 0, Y is a singular (or smooth) solution, and so can be also 
V. Moreover, in the case that N is transversal to a characteristic smooth vector 
field C : Ek — > TEk, then V = Uag]-e £ [ ( / ) a(^) is a smooth solution of Ek passing 
through N. Under suitable conditions on the Cauchy integral manifold, the vector 
field used to build the solution can be an infinitesimal symmetry. (For a detailed 
proof see the one of Theorem 2.10 in [55].) 

Finally, whether Ek — > E^—i = ~Kk,k-i{Ek), is an affine subbundle of iTk,k-i ■ 
J„(W) — > J k ^ 1 (W), with associated vector bundle 7ffe_i : gk — > E^-i, then also Ek 
is a strong retract of J*(W), so we can reproduce above strategy used to build a 
solution passing for N, without the necessity to prolong Ek, (if it is not differently 
required by the structure of this equation). □ 

Example 3.21 (Fourier's heat equation). Let us consider the second order PDE 

(32) (F) C JD 2 {W) C Jf (IF) : F = u t - u xx = 

on the fiber bundle n : IF = M 3 — > R 2 , (t,x,u) >— > (t, x). In this case (F) — > 
TV2,i(F) = (-F)-i = JD(W) is an affine fiber subbundle of the affine fiber bundle 
7T2,i : JD 2 (W) — > JD(W), with associated vector bundle identified with the symbol 
<72' C = Qttduu + Qxdutx G <72- So we can apply Lemma 3.20(3). Therefore, if N C 
(F) is a (compact) 1-dimensional integral manifold, diffeomorphic to its projection 
t^2,o(N) = X C IF, we can find solutions V C (F), passing from N. In particular, 
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whether N is diffeomorphic to a smooth space-like curve u = h{x), (at t = 0), we 
get that N is the image of a mapping, say 7 : / = [0, 1] C R — > (F), iff 7*02 = 0, 
where O2 =< dF, u> = du — u t dt — u x dx, u>t = dut — u u dt — u tx dx, uj x = du x — u x tdt — 
u xx dx > is the Pfaffian contact ideal encoding solutions of (F). Then one can see 
that the integral curve 7 is represented in coordinates (t, x, u, u t ,u x , Utt, Ut x , U xx ) on 
JD (W), by the following equations: 
(33) 

J £07 = 0; .to 7 = 2;; uo-/ = h(x); u 4 o 7 = 0(ir); u x oj=^(x); 
\ u xx ° 7 = jj^( x )'i u tx o 7 = ^t(x); u«o7 = k(x). 

where k(x) is an arbitrary smooth function. Then, by considering that the Fourier's 
heat equation is a formally integrable and completely integrable PDE, we can see, 
by taking the first and second prolongations of (F), that must necessarily be k{x) = 
(x) . In fact, from the first prolongation we get that must be Utt = u xxt . From 
the second prolongation of (F) one has u tX x = u XX xx- Since must be u xx t = Ut 
we get u tt = u XX xx- By restriction ont = 0. one has u tt \t=o = u xxxx \ t =o = h. 
So we can build solution also by using Lemma 3.20(1). For example, if we are 
interested to a solution Vs C (F ), obtained by means of a rigid propagation of the 
initial Cauchy space-like integral curve (33), one can see that must necessarily be 
h(x) = ax + (3, with a, (3 £ M., x G [0, 1] C R. In other words such type of steady- 
state solution, u(t, x) = ax + ft, determines also the admissible structure of the 
integral initial Cauchy line T C (F) t= o. In such a case this must be given by the 
following parametric equations in JD 2 (W): 
(34) 

407 = 0; 107 = 1; u o 7 = ax + /3; ut ° 7 = 0; u x o 7 = a; 

u xx °7 = 0; M(i°7 = 0; utt °7 = o. 

Remark that C = dt is the smooth vector field propagating the initial Cauchy curve 
r in Vs , i.e., £ is the characteristic vector field ofVs- This is not a characteristic 
vector field for the Fourier's heat equation and neither it is a characteristic vector 
field for the steady-state sub-equation {u xx = 0}. Instead, £ — dt is an infinitesimal 
symmetry for (F). (See Lemma 3.20(2).) Let us also emphasize that the above 
steady- state solution is not the unique solution passing for V . In order to see this 
it is enough to generalize the concept of solution and to consider also weak-singular 
solutions. In fact, let us find regular perturbations of the steady-state solution Vs 
that at t = 0, in correspondence of the boundary points {A, B} = dT of T, have 
values respectively a and a + b. These can be obtained by considering the Fourier's 
heat equation, that is a linear equation, as its linearized at the steady-state solution 
Vs, (F)[Vs] = (F), i.e., the equation for perturbations 0/V5. Then we get: 

h b — 6cos(/i) 



(T) 



(35) v{t, x\n, a, b) = e~» t 



sm(^j,x) + b cos(fj,x) 



sin(/i) 

v(t,x\fj.,a,b) depends on an arbitrary positive parameter (i > 0, and two other 
parameters a, b E R, and has the following limits: 

lim v(t, x\/i, a, b) = 0; lim v(t, x\[x, a, b) = ax + 6; Iim v(t,x\fx,a,b)=0. 

t— y+00 fi— >0 a,b— >0 

This means that with respect to these perturbations the steady-state solution is 
asymptotically stable, and all such perturbations produce deformations of the steady- 
state solution, with deformation parameters just (//, a, b). Let us denote by V[jJL, a, b] 
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the 2- dimensional integral manifold, contained in (F), representing the deformed 
steady-state solution by v(t,x\fj,,a,b), < x < 1. More precisely V[fi, a, b] is the 
integral manifold of (F) corresponding to the solution u(t,x\(J,,a,b) = ax + /3 + 
v(t,x\fj,,a,b). Set V = Vs[jV[fi,a,b]. V is a weak-singular solution of (F), that 
for t = passes for S^ 2 -* = (J a, b]( 2 \ a 1-dimensional admissible integral 
manifold of (F), such that n2,o{'£ ) = S = r (J a, b], where T[/j,,a,b] is the 
curve identified by V[/i, a, b] for t = 0. 3 is a space-like curve contained in W , that 
is projected, by means of it, on the interval [0,1] of the x-axis. The 2- dimensional 
manifold Y = ir 2 ,o{V) = 7r 2 ,o(Ks) U ^2,o(V[fj., °> °}) = Y S U Y \M> a -> b ]> passes for S. 
Furthermore V converges, for t —> +oo, to the steady-state solution Vs. 19 
Let us, now, show also an explicit construction of solution for (F) by means of 
the method of the retraction given in Lemma 3.20. For example, let f(P) = V C 



Ji(W), P c 



be a 2- dimensional integral manifold representing the smooth 



function u(t,x) = ^ /n >ot n a n (x), passing for the integral curve given in (33). The 
coefficients a n (x) are suitable functions on x £ [0, 1], such that the series converges 
in t € M. + , and such that ao(x) — h(x), a\{x) = h xx , a^ix) = \h xxxx . These last 
conditions assure that V passes for the integral curve (33). The integral manifold V 
is not a solution of (RF) since there are not restrictions on the the other coefficients 
a nt n > 3. The parametric equations of V are given in (36(a)). 
(36) 

t°f = t 

xof = x 

« / = Eo<n<2 + En>3 «"«»(*) 

«* / = Ei<„< 2 ^J^W + E„> 3 nt-^a^x) 
«. o / = E <n<2 viSW + E„> 3 t^(x) 
«« ° / = Eo<n< 2 SSS^W + E„> 3 

° / = E!<„<2 i- 1 ^ (*) + E„> 3 „ ^-^W 



■V 

J/ 



(b)< 



t o r o f = t 

x o r o f = x 

u o r o / = u(t, x) 

u t o r o f = u t (t, x) 

u x or o f = u x (t,x) 

u tt o r o / = u a (t,x) 

u tx ° r o / = S fx (t,a;) 



i3 n(n - l)f 2 a„(x). 

Taking into account that on (F) we can consider the following coordinate functions 
{t,x,u,Ut,u x ,Uu,Utx}> it follows that the retraction mapping r(f(P)) = r(V) 
V c (F) has the parametric equation given in (36(b)), where the function u(t,x) is 
determined starting from the function u{t,x), by imposing the condition to belong 
to (F). So we get u(t,x) = E n > «n(^), with a n = ^^(x) 20 V C (F), 



19 Singular solutions, like those described in this example, are very important in many physical 
applications too, since they represent complex phenomena related to perturbations of some fixed 
dynamic background. For example, for suitable values of the parameters a and b in (35), manifolds 
Y[fi, a, b] intersect Y$ along common characteristic lines, and the singular solutions V are piecewise 
Z2-manifolds. (For complementary informations on such singular manifolds, and singular solutions 
of PDE's, see also [13, 42, 55, 62].) 

This function u(t, x), is well denned and limited in all (t, x) £ M. + X [0, 1], since the function 
h(x) is smooth in [0,1] and with ||^#0)| < C £ 1, Vi 6 [0,1]. In fact one has £„> i™S„ < 



C^2 n>0 b n - The last series converges, with convergence radius r = 00 since 



lim„_ 



lim n _ 



t 

n + l 







Therefore the series 5Zn>o*" an ( x ) * s convergent as it 



is absolutely convergent. 
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represents the solution u(t,x). This can be considered as a deformation of u(t,x). 
In fact, the functions defined in (37), is just the explicit deformation connecting 
u(t, x) with u(t, x). 

u(t,x\X) = 52n>o tn ®n(M x )i a n (X\x) = a n {x) - X[a n {x) + a n {x)} 
(37) | u(t,x\0) = j: n > t n a n (x) 

{ u(t,x\l) = E„> 

The integral manifolds V\ C J%(W), A G [0,1], corresponding to u(t,x\X), are 
not contained in (F) for all A G [0, 1], hence are not solutions of (F). Therefore, 
V C (F) is a solution for the corresponding Cauchy-problem, obtained by means of 
the retraction method given in Lemma 3.20. The parametric equation ofV is given 
in (38) 



( torof = t 



(38) 



x o r o f = 


X 


u o r o f = 


4-n 1 d 2n h („\ 
Z^n>0 L n \ dx 2n W 


u t orof- 


- Z^n>l 1 (n-1)! dx 2 " W 


MjOro/: 


_ ,nl d 2n + 1 h 1 \ 
~ l^n>0 L n! dx 2 ^ 1 W 


u t t°rof 


— fn-2 1 d 2 "h 1 \ 

- l^n>2 1 (n-2)! dx 2 " W 


u tx or o f 


_V f 1 - 1 1 d 2n+1 h( \ 
l^n>l L (n-1)! dx 2 ™+ 1 W 



> x e [0, 1] 



Let us emphasize, that also in this case the solution so obtained is not unique. 
In fact, similarly to the previous case, where we have considered weak-singular 
solutions by means of perturbations of Cauchy data for the steady state solution 
Vs, now we have the following weak-singular solution V = V [jV[(i, a,b] C {F), 
where V[/i, a, b] is the 2-dimensional integral manifold identified by the solution of 
(F) given in (39). 
(39) 

1 d 2n h , s _„2 t \a + b - 6cos(/i) 



l(t,x\n,a,b) = 



n>0 



i! dx 2 '- 



-(x) + e 



— u t 



sin(/j) 



siii(/ia;) + b cos(^ix) 



Let us denote by T[/j,,a, b] C Wt=o the space-like curve identified by u(0,x\fi,a,b). 
Then V[n, a, b] passes for 

rW[n,a,b} C (F)t=o, 
and the singular solution V passes for the space-like curve 
5(2) = f( 2 )|Jf (2 W,6] cVc(F), 

where 

E = f\Jf[fx,a,b] C Wt=o 
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is the space-like curve containing the fixed Cauchy data, i.e., the curve T. Further- 
more the weak-singular solution V asymptotically converges (t — > +00) to the regu- 
lar solution V. Since the singular solution V depends on the parameters (fi,a,b) G 
R + x R 2 , the Cauchy problem identified by the curve T has more than one solution. 

Let us, now, Nq,Ni C E k be two closed compact [n — l)-dimensional admissible 
integral manifolds of Ek- Then there exists a weak, (resp. singular, resp. smooth) 

solution V C E k , such that dV = N (j Ni, iff X = N (j Ni G [0] G Sl%l ltW , (resp. 
X G [0] G fifii s , resp. X G [0] G &n-i)- On tnc other hand there exists such a 
smooth solution iff X G [0] G f2„_i and X has zero all its integral characteristic 
numbers, i.e., are zero on X all the conservation laws of Ek- Since these last can be 
identified with the Hopf algebra H n -i(Ek) = H n -i(E 00 ), where A^ is the infinity 
prolongation of Ek, it follows that the quotient H n -i(E 00 ) /M.^ 1 ™- 1 measures the 
amount of how the conservation laws of Ek differ from the crystal conservation 
laws, identified with the elements of the Hopf algebra K^™- 1 . □ 

Definition 3.22. We define crystal obstruction of Ek the above quotient of alge- 
bras, and put: cry(Ek) = H n -i(E 00 )/R n ' 1 - 1 . We call 0-crystal PDE one Ek C 
J*(W) such that cry(E k ) = 0. 21 

Corollary 3.23. Let E k C J*(W) be a 0-crystal PDE. Let N ,Ni C E k be two 
closed compact (n — 1)- dimensional admissible integral manifolds of Ek such that 

X = Nq (J Ni G [0] G Qn-i- Then there exists a smooth solution V C Ek such that 
dV = X. 

Example 3.24 (The Ricci-flow equation). The Ricci-flow equation 
(40) F tJ = (dt. 9lJ )- K R l3 =0 

on a Riemannian n-dimensional manifold (M,g), can be encoded by means of a 
second order differential equation (RF) C JD 2 (E) C J%{E) over the following 

fiber bundle tt : E = K x S%M — S- R x M, where S%M C S%M is the open 
subbundle of non- degenerate Riemannian metrics on M . In [60] we have calcu- 
lated the integral bordism group of the equation (RF) C J|(£?). In particular if 
M is a 3- dimensional closed compact smooth simply connected manifold, we get 
^Is^uf' — ^3^^ — ^2- Thus (RF) is not an extended 0-crystal PDE. Taking 
into account exact commutative diagram (27), we get also the short exact sequence 
— > A3 ' — > fig — > %2 — > 0. Taking into account Example 2.19 we can consider 

f2g as an extension of a subgroup of the crystallographic group G(3) = Z 3 x Z2 
or G(2) = Z 2 x %2- Therefore the integral bordism group of the Ricci-flow equa- 
tion on S 3 is an extended crystal PDE, with crystal group G(2) = Z 2 x Z2 = p2 
and crystal dimension 2. Furthermore, from the exact commutative diagram (27) 

(RF) (RF) 

we get also the following short exact sequence 

Taking into account that O3 = 0, we can have cry(RF) ^= 0. On the other hand let 
us consider admissible only space-like integral Cauchy manifolds satisfying the fol- 
lowing conditions: (i) They are diffeomorphic to S 3 or to M , assumed any smooth 
3- dimensional Riemannian, compact, closed, orientable, simply connected manifold; 

21 An extended 0-crystal PDE E k C J*(W) does not necessitate to be a 0-crystal PDE. In fact 
E k is an extended 0-crystal PDE if fif?^ w = 0. This does not necessarily implies that £l I ^ k _ 1 = 0. 
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(41) 



► , 1 < i,j,h,k < 3. 



(ii) M is homotopy equivalent to S 3 . Such integral manifolds surely exist, since we 
can embed in E both manifolds M and S 3 and identify these, for example, with 
space-like smooth admissible integral manifolds of the subequation (RF) t C (RF), 
where (RF) t = tt^ Wj t £ M.. Here 7T2 is the canonical projection (RF) — > K, 
induced by n 2 : JD 2 {E) — > E x M , i.e., n 2 = pr\ o 7r 2 : JD 2 {E) — > R, where 
pr\ is the canonical projection pr\ : M X M — > M. In /act, in (RF), (dt.gij) is 
solved with respect to Rij . More precisely, starting from a 3- dimensional compact 
closed, orientable, simply connected Riemannian manifold (M, 7), we can identify 
a space-like integral Cauchy 3- dimensional manifold No C (RF)t=t =o, diffeomor- 

phic to its projection tt2,o(Nq) = Yq C Wt=t =a — S®M , by means of the mapping 
f : M —> JD 2 (E) defined by the parameter equation in (41). 

i o / = t = 
x k o f = x k 

9i] ° f = la 

9ij,t of = K>Rij(l) 
9ij,h f = HjJi 
9ij,th ° f = KRij(j),h 
9ij,hk ° f — Jij,hk 

9ij,tt °f = 

where {t, x k , g i3 , g ijit , g ijih , gij,th, gij,hk, 9ij,tt}i<i,j,h,k<3 are local coordinates on JD 2 (E), 
<&ij{pi), i,j £ {,2,3}, are known analytic functions of 7^ and its derivatives up to 
fourth order, symmetric in the indexes. In fact, since (RF) is a formally integrable 
and completely integrable PDE [60] , from its first prolongation we get 

{9ij,tt = K Rij(g),t,gij,th = K-Rij {g),h}l<i,j,h<3, 

where Rij(g)j and Rij(g)^ denote the first prolongation of Rij (g), with respect 
the t-variable and x h -variable respectively. Since Rij(g) = Rij(g r s, grs,h, grs.hk), 
i,j £ {,2,3}, are analytic functions, we get also the following analytic functions 

Rij(9),t — Kij{9rs, 9rs,h, grs,hk, grs,t, grsM, grsMkt), 1 < i, j < 3 

where Kij are analytic functions of their arguments. Taking into account the expres- 
sion of (RF) and second order prolongation, we get the initial values (^-dimensional 
Cauchy integral manifold), given in (42). 



(42) 



grs,h\t=Q — IrsJi 
grs,t\t=0 = nRrs(l) 
9rs,th\t=0 — KRrs,h\t=0 
gra,hk\t=0 = lrs,hk 
grs,tt\t=0 = ®ij(l) 



Kr,sji,k<3 



where R rs ,h is the first prolongation of R rs with respect to the coordinates x h . Tak- 
ing into account that one has the following functional dependence R rs = R rs (gij, gij, p , gij, P q) 
we get 

Rrs,h = (dg ij .Rrs)9ij,h + (dg lJ ' p .R rs ) gi ^ ph + (dg %3 > pq .R rs )g %] , pqh . 

As by-product we get that R rs ,h\t=a are functions of 7^ and their derivative up 
to third order, that we shortly denote by R rs ,h\t=o- Furthermore, from the first 
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(43) 



prolongation of (RF) we get also 

g rs ,tt = nR rs ,t = K[(dg' ] .R rs )gij, t + (9g' :hp .R rs )gij, P t + (dg lhpq .R ra )gij, pq i\ 
= K[(dgv .R T3 )K,R tJ + {dgW.Rrs^Rij* + (dgV™ .Rrs)nRi3, Pt } 

Then taking t = 0, we get that g r s,tt\t=o functions that depend on jij and their 
derivatives up to fourth order, that we shortly denote by $ij (7) . 



l<r,s,h,k<3 



(44) 



JD 2 (E) ^-^ JD 1 (E) ■ 



- S°(U x M) ® S°(M) *- vTJD 2 (E) w^vTJD^E) 



■ vT(RF) ^ 7rJ lt )T( J RF)_i 



' (RF) ■ 



■ (RF)- 





Similarly one can identify the Riemannian manifold (S 3 ,^/), with another space- 
like integral Cauchy 3- dimensional manifold Ni c (RF) t= t 1 ^o, diffeomorphic to 

its projection ^2,0(^1) = ii C Wt = t 1 = S®M. Let us emphasize that, fixed the 
space-like fiber (RF) t , above integral Cauchy manifolds are uniquely identified by the 
Riemannian manifolds (M, 7) and (S 3 ,j), respectively. For such integral manifolds 
Nq and Ni, necessarily pass solutions of {RF), (hence they are admissible). 22 We 



In general a steady state solution of (RF) is not admitted since this should imply that 
(M, 7) is Ricci flat. Furthermore, regular solutions gij(t,x k ) = h(t)gij(x k ), with separated time 
variable from space ones, imply that (M, 7) is an Einstein manifold. In fact Rij (g) = Rij (hg) = 
Rij(g). The Ricci-flow equation becomes ht(t)g~ij(x k ) = nRij(g). Therefore, must be ht(t) = 
ui = KRij(g)/g~ij(x k ), with ui £ R. By imposing the initial condition gij (0, x k ) = r yij(x k ), we 
get that must be h(t) — uit + 1 and g~ij(x k ) = ■yij(x k ), hence Rij(-y) = -fijUi/K. Vice versa, if 
the Ricci flow equation is considered only for Einstein manifolds, then solutions with separated 
variables like above, are admitted. This means that in general, i.e., starting with any (M, 7), 
we cannot assume solutions gij(t,x k ) with the above separated variables structure, even if these 
solutions "arrive" to S 3 , that is just an Einstein manifold. The same results can be obtained by 
considering metrics g(t, x k ), obtained deforming r y(x k ), under a space-time flow <f>x of M X R, i.e., 
{t o if>\ = i + A, x k o 4>\ = <j>\ (t, £ l )}i<i,fc<3, and with initial condition {0q = x k }k=i,2,3 ■ I n f ac t 
if we assume that (f> r x = h(X)<j> r (x k ) , then the Ricci flow equation becomes as given in (45). 

u<\ fr 2 M _ i: VM R rsh)°<t>x) _ , . r p+^ f (a) h(X) - ±^h(\) = \ 
y ' h\\) ##(>.(7Wa) l( b ) ^[M7)-;>.]^ = J 

The integration of the equation (45(a)) gives h(X) = Ce ±x ^" . Furthermore, if (M, 7) is an 
Einstein manifold, i.e., there exists fi S M, such that i? r s(7) = M7rs, then taking u) = fiK, 
one has the following solution, uniquely identified by the initial condition, (up to rigid flows): 
<f>x = e ±A v / A" ? £ r . If (M, 7) is not Einstein, or equivalently, assuming ui ^ Kfi, we see that the 
solutions of equation (45(b)) are <j> r = a r 6 R, since the metric 7 rs = R rs (-f) — ^7rs is not 
degenerate, hence must necessarily be cf> r = 0. But such a flow does not satisfy initial condition 
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shall prove this by considering Lemma 3.20(2). First let us note that, similarly to 
the Fourier's heat equation, the Ricci-flow equation is an affine subbundle of the 
affine bundle W2,i ■ JD 2 (E) — ¥ JD 1 (E), since the fiber 7r^"J(g) = (RF) q , for q G 
(RF)-i, is an affine space. In fact, T q (RF) q = (92)13, and the vector space (92)9 is 
constant on the fiber q S (RF)—i. The situation is shown by the exact commutative 
diagram in (44) and by the fact that a vector field £ = Ct/^dap belonging to the 
symbol must satisfy equation (46). (One has used the usual space-time indexes 
numbering for the coordinates (t,x k )i<k<3 = (a; Q )o<a<3 on M x M.) That equation 
is constant on the points of the fiber (RF) q , since it depends only on the derivatives 
of first order gij.h and zero order g^. 



(46) Cf(dg r a s p .F ij )(q) = 0,qG(RF). 

This follows soon from the expression of the Ricci tensor as a differential polyno- 
mial. See equation (47). 



(4-7*) R — n r PR / Rrjnp — 2^ r P^n + Qjn.rp 9rn.jp Qjp.rn) \ 

{ ' 9 nnp \ +g t '(\jn,s\[rp,t]-\jp,s\[rn.t\) } 



[ij, k] = \{gikj + gjk., - g,j.k) 
g rp = MM 

\g\ = det(g rp ) 

[g r p] '■ algebraic complement of g rp . 



Since the derivatives of second order gij^hk appear in (47) multiplied by zero-order 
terms only, it follows that equation (46) is constant on each fiber (RF) q . As a 
by-product it follows that we can apply the retraction method given in Lemma 3.20 
directly on (RF). 

With this respect, let us see more explicitly, and less formally, how we can build, 
a smooth solution of the Cauchy problem (RF) identified by the 3- dimensional in- 
tegral manifold given in (42). Since, fixed the Riemannian manifold (M, 7), we 
uniquely identify a 3-dimensional integral manifold Nq C (RF), diffeomorphic to 
M, (resp. S 3 ), then in order to build a solution passing for Nq, (resp. N\), it 
is enough to prove that we are able to identify a time-like vector field, tangent to 
(RF), transversal to Nq that besides the tangent space TNq generate an integral 
planes of (RF). In fact, integral planes Lq C (E2) q = L q Q)(g2) q C T q (RF), 
q G (RF) + 4, q = 7^2(9), are generated by the following horizontal vectors Ca(<l) = 
[dx a + Y.o<\0\<29ij,apdg' i: i< l3 ]- q( z> RF ) +i € T q (RF) . (\/3\ denotes the multiindex length 
and < a < 3.) In (48) is given a more explicit expression of the horizontal vectors 

Co- 
rn 

Co = dt + gutdgV + gij.thdg^- 11 + ga.ttdg^ + gi j,thkdg i3 ' hk + gij.tthdg^ + g^utdg^ 1 
gij.t = nRij{"f) , gij^h = Ki?ij,h(7) 

Sy.a = >tRij,t(i) = « Idg" .%)reii rs (7) + (dg"-" .R l:i )KR rs<p (y) + {dg TS ™ .i? IJ )K.i? rs , M ( 7 )] 

gij.thk = K-Rij,hk(l) ; gij,tth = KRij,th{"f) J 9ij,ttt = K Rij,tt{l) 



{o5q = x r }r=i,2,3- In conclusion a metric gij(t, x ), obtained deforming 7 with a space-time flow 
<f>\, where 4>\ = h(\)<j> r (x k ), is a solution of the Ricci flow equation iff (M,7) is Einstein, or 
Ricci-flat. This last case corresponds to take ui = in equation (45) and has as solution the 
unique flow, up to rigid ones, <p r x = x r . 
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Figure 1. Solution Cauchy problem for (RF), represented by the 
integral manifold V, envelopment manifold, generated by local so- 
lutions represented by manifolds V q tangent to TVo, Vq £ No, and 
identified by means of 4-dimensional integral planes L q , for any 

q £ A^q 1 -*. V ou tside C JD 2 (E) represents a 4-dimensional integral 
manifold, passing for N , but V outS ide <t (RF). Voutside can be de- 
formed in V, taking fixed Nq. (RF) and (RF) + i are represented 
with (yellow) frames, containing respectively Nq and Nq . 



ForqE (^Vo) (4) = D 6 -f(M) C [(RF) +4 } t = C (RF) +4 , at the corresponding point 
q £ N , one has ( = Qs) + C(T) G ( E 2)g. ««ift C(s) = X k ( k (q) £ T 9 iV , and 
C(t) = -^°Co(<?) transverse to Nq. Thus the 4-dimensional integral manifold V q , 
tangent to such an integral plane, is tangent to Nq too, and has the vector £o as 
characteristic time-like vector at q £ Nq. By varying q in Nq we generate a 4- 
dimensional manifold V that is the envelopment manifold of the family {V q } q ^N 
of solutions of (RF), formally defined asV = UgeJV -^-o.q, where Xq^ is the integral 
line transversal to Nq, starting from q, tangent to Co(l) and future- directed. Thus V 
is a line bundle over Nq, containing Nq, with T q Nq C T q V , Vq £ Nq. This is enough 
to claim that V is contained in (RF) and it is a time-like integral manifold of (RF), 
whose tangent space is an horizontal one. In other words, V is not only a viscosity 
solution, but just a solution of the Cauchy problem identified by Nq (resp. N±). 
(See Fig. 1 where the 3- dimensional manifold Nq is reduced to dimension 2 (figure 
in the left- side) or dimension 1 (figure in the right-side) for graphic necessities.) 
Let us emphasize that, on Nq, (resp. N\), the components of the characteristic 
vector field C,q in (48) are uniquely characterized by the metric jij and its deriva- 
tives up to sixth order. This proves that in order to characterize a solution of the 
Cauchy problem, identified by Nq C (RF), it is necessary to consider the prolon- 
gation (RF) +4 , of (RF), up to fourth order. Then the 4-dimensional manifold 
V = Ute[o e[ 'Pt(No) C (RF), obtained by means of the local flow <f>t generated by 
Co = d<p, is a local integral manifold, solution of the Cauchy problem identified by 
Nq. (So Co is the characteristic vector field for such a solution.) 24 Let us remark 

^Generalized solutions of PDE's, called viscosity solutions were introduced by Pierre-Louis 
Lions and Michael Crandall in the paper [36]. Such solutions do not necessitate to be solutions, 
but are envelope manifolds of solutions. 

24 Let us emphasize also that an integral 4-plane, where the components gij,g, < < 3, 
satisfy conditions in (42), but not all conditions in (48), has the corresponding integral manifold, 
say V, that passes for Nq, but it is not contained in (RF). (This is represented by V ou tside m 
Fig. 1.) Let us denote the corresponding metric with g. Therefore the retraction method imposes 
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(50) 



that such a solution does not necessitate to be smooth. In fact in this process we 
have used only the prolongation of (RF) up to fourth-order, i.e., we have considered 
only the projections given in (49). 

(49) (RF) +4 (RF) +3 (RF) +2 ^ (RF) +1 (RF) . 

However, we can obtain V as a continuous manifold. In fact, let gij{x a ) denote 
a local solution identified by Lq- = T q V q C T q (RF), and T q No C Then the 

time-like integral line -Xo,g C V q , tangent to Co, is represented by a local curve 
£- q ■ [0,e g [c M — > (RF), given by the parametric equation (50) in JD 2 (E). 

to£ = t 
x k o £ = x k (q) 
9ij °£ = 9i](t,x k (q)) 
9ij,t°€ = KR l:j {g(t,x k ))\ x k =x k {q) 
9ij,h °£ = 9ij,h(t>% )\x k =x k (q) 
9ij,hk ° £ = 9ij,hk(t> ^ )lx fe =s fe (g) 
9ij,th °£ = K Rij,h(g{t,x k ))\ x k =x k( q } 

9ij,tt°t = ^ 2 ilo<\ p \< 2 \( d 9 rs,p -Ry)Rrs, P {g(t,x k ))]\ xk=xk{q] 

Since Co continuously changes on Nq, it follows that curves Xq^, q G Nq, continu- 
ously change with q G Nq, if the solutions gij(t,x k ) continuously change with q G Nq 
too. This is surely realized by the fact that points q' G iVo, near to q, come from 
points q 1 G Nq , near to q, if ^3,2(9) = Q- In f a ct, Nq = Nq and ^3,2 (iVg ) = No, 
thus 7T3 : 2 1 Tv-t 1 ) : -^o 1 ' — ^ -^0 i s necessarily a continuous mapping. Moreover, we can 

continuously transform any (local) solution gij(t,x k ) of {RF) into other ones by 
means of space-like (local) 1-parameter group of diffeomorphisms (f>\, A G [0,e[c K 
of M . In fact, such transformations are locally represented by the following func- 
tions {to(f)\ = t,x k o(f)\ = <^(x 1 )}. As a by-product it follows that (RF) is invariant 
under such transformations. In fact (4>*\9)ij,t~ K Rij(4 l \9) = fixlilijj — KRij(g)] = 0. 
Thus we can continuously transform an integral manifold V q , into other solutions 
along the space-like coordinate lines of Nq, passing for q and identify, in the points 
q' G Nq, the time-like curves Xq^ that result continuously transformed of XQ q . In 
this way V = UggAf ^o,q is a smooth manifold in a suitable tubular neighborhood of 
Nq x [0,e[, that is the integral manifold of a metric gij(t,x k ), of class C 3 , solution 
of the Ricci flow equation (RF) . This proves that the envelopment manifold V is 
more regular than a viscosity solution. 

In order to obtain a smooth solution of the Cauchy problem given by the integral 
manifold Nq it is necessary to repeat the above process on the infinity prolongation 
(RF) +oc C JD°°(E). In fact, also on (RF) +!X1 , the ^-dimensional integral manifold 
Nq = D°°^(M) = Nq 00 ^ G (RF) +oa is uniquely identified by (M, 7), and identi- 
fies also a unique transversal time-like characteristic vector field Co > tangent to 



also to gij-,/3, with < \f)\ < 3, to satisfy conditions reported in (48). The relation between 
V = V ou tside C JD 2 (E), and V C (RF), can be realized with a deformation, g\, connecting the 
corresponding metrics g and g. More precisely, g\ = g-\-X[g — g\, A £ [0, 1]. The integral manifolds 
V\ C JD 2 (E), generated by g^, are not contained in (RF) for any A £ [0, 1], but all pass for Nq. In 
fact, since = 9ij,p+tyij,p-9ij,p], I/ 3 ! > 0, we get that (g\)ij^\ No = gij,p\N = 9ij,p\N , 

for < \P\ < 2, but (g\)ij : p |jv do not satisfy conditions in (48) for \f3\ = 3. This is just the 
meaning of the retraction method considered in the proof of Lemma 3.20. 
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(RF) +oa . More precisely, Q = Co + J2\ a \>2 9ij,tadg v a , with g ij>ta = nR ijta (j), 
where a = 0, 1, 2, 3. 

The situation is resumed in the commutative diagram (51). 

(51) T(RF) +X ~ ^TN^ ] 



T(RF) +X ^—^(RF) +00 * 3 [(iJF) +oo ] t=0 < jjyWc JD°°(S°M) 




HxM -« jjlf = M = A/ 



Similarly to what made in the Fourier's heat equation, (see Example 3.21), we 
can prove that to the above solution one can associate weak-singular ones by using 
perturbations of the initial Cauchy data. In fact, if V C (RF) is a regular solution 
passing for the integral manifold No, we can consider perturbations like solutions 
v : R x M ->• E, of the linearized equation (RF)[V] C JD 2 (E). Since also (RF)[V] 
is formally integrable and completely integrable, in neighborhoods of No there exist 
perturbations. These deform V (background solution) giving some new solutions 
V C (RF). Then V = V(J V c (RF) is a weak-singular solutions of the type just 
considered in Example 3.21 for the Fourier's heat equation. 2 ^ 

Now, for any of two of such integral manifolds, No and N\ , we can find a smooth 
solution V bording them, V — No (J N\ iff their integral characteristic numbers are 
equal, i.e. all the conservation laws of (RF) valued on them give equal numbers. By 
the way, under our assumptions we can consider No and Ni homotopy equivalent. 
Let f : N\ — > No be such an homotopy equivalence. Let uj be any conservation law 
for (RF). Then one has equation (52). 



(52) < [N ], M >= f u = [ f*oj =< [N{\, lf*Lo] >=< [N^, [to] >= [ u. 

JNq JNi JNx 

So any possible integral characteristic number of No must coincide with ones with 
Ni and vice versa. Thus we can say, that with this meaning of admissibility (full 

2 ^This generalizes a previous result by Hamilton [27] , and after separately by De Turk [19] and 
Chow & Knopp [18], that proved existence and uniqueness of non-singular solution for Cauchy 
problem in some Ricci flow equation. Let us also emphasize that our approach to find solutions 
for Cauchy problems, works also when JVq C {RF) is diffcomorphic to a 3-dimensional space-like 

submanifold of W , that is not necessarily representable by a section of Et=o — S®M — > M. 
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admissibility hypothesis^ on the Cauchy integral manifolds, one has cry(RF) = 0, 
i.e., (RF) becomes aO-crystal. Therefore there are not obstructions on the existence 

of smooth solutions V of (RF) bording No and Ni, dV = -/Vo(JiVi, i.e., solutions 
without singular points. 

This has as a by-product that M and S 3 are homeomorphic manifolds. Hence 
the Poincare conjecture is proved. With this respect we can say that this proof 
of the Poincare conjecture is related to the fact that under suitable conditions of 
admissibility for the Cauchy integral manifolds, the Ricci-flow equation becomes a 
0-crystal PDE. 26 

Example 3.25 (The d'Alembert equation d dx °Qj = on the 2-dimcnsional torus). 
The d'Alembert equation on a 2- dimensional manifold M can be encoded by a 
second-order differential equation 

(53) (d'A) 2 C JD 2 (W) : {uu xy - u x u y = 0} 

with W = M x R. In [60] we have calculated the integral bordism groups of such 

an equation. In particular for M = T 2 , the 2-dimensional torus, one has = 

2^2 ©^2- Taking into account Example 2.20 we see that is isomorphic to 

the crystal group G(2) = T? x D± = p4m. Therefore, (d' A) 2 on T 2 is an extended 
crystal PDE, with crystal dimension 2. 

Example 3.26 (The Tricomi equation on 2-dimensional manifolds). In [60] we 
have considered the integral bordism groups of the Tricomi equation (T) : 

(54) (T) C JD 2 (W) : {u yy - yu xx = 0} 

defined on a 2-dimensional manifold M , i.e. with W = M x R. For example, on the 
2-dimensional torus T 2 one has = Z2 © Z2. Furthermore on M.P 2 we obtains 
= Z2. Thus the Tricomi equation on T 2 , (resp. S 2 ), is an extended crystal 
PDE with crystal group p4m, (resp. p2), and crystal dimension 2, (resp. 2). 

Example 3.27 (The Navier-Stokes PDE). The non isothermal Navier-Stokes equa- 
tion can be encoded by a second order PDE (NS) C JD 2 (W) on a 9- dimensional 
affine fiber bundle n : W —> M on the Galilean space-time M . In Tab. 1 is reported 
its polynomial differential structure. 

There {x a ,x k ,p, 9} are fibered coordinates on W, adapted to the inertial frame, and 
G\j are the canonical connection symbols on M . Furthermore A C IR^x 1 , x 2 , x 3 ]] is 
the algebra of real valued analytic functions of (x k ). 27 We have proved in Refs.[55, 

2 ®Thc proof of the Poincare conjecture given here refers to the Ricci-flow equation, according 
to some ideas pioneered by Hamilton [26, 27, 28, 29, 30], and followed also by Perelman [43, 
44], However the arguments used here are completely different from ones used by Hamilton and 
Perelman. (For general informations on the relations between Poincare conjecture and Ricci- 
flow equation, see, e.g., Refs.[5, 16, 17] and papers quoted there.) Here we used our general 
PDE's algebraic-topological theory, previously developed in some works. Compare also with our 
previous proof given in [2, 3], where, instead was not yet introduced the relation between PDE's 
and crystallographic groups. 

ffijfx 1 , . . . , x n }} denotes the algebra of formal series JT, j in o-i 1 ...i rl (a: 1 ) 11 ■ ■ ■ (1")'" , with 
ai 1 ...i n S IR. Real analytic functions in the indeterminates (a; , . . . , x n ), are identified with above 
formal series having non-zero converging radius. Thus real analytic functions belong to a subal- 
gebra of R^rr 1 , . . . , x n \\. This last can be also called the algebra of real formal analytic functions 
in the indeterminates (a; 1 , . . . , x n ). 
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Table 1 . Completely integrable Navier-Stokes equation: (NS) C 
JD 2 (W) denned by differential polynomials. 



(A) 


: F° = ± k G] k + ±16? = {F° aA[± k ,±i\) 

(continuity equation) 


(B) 


: F° = x k (dx a .G] k ) + x k a G] k + x\Jf = (^eA^,*- J) 

(first prolonged continuity equation) 


(C) 


: Fi = x s Ri + x s x l pGi + x s x j s p + px J + x k S J k + x\ s T ls 

+Pig ij + p(dx t .f)g lJ = (F j eA[x",xi,xi 3 , P i]) 

(motion equation) 


(D) 


: f 4 = e oP c p + pc p x k e k + e^ s + x k ±pw kp + x k ± s p w p ks + >vr = o 

(F i £A[x k ,x s p ,0o,Ok,8u]) (energy equation) 


Functions belonging to ACRffa; 1 ,x 2 ,x 3 ]] 


Ri= 


xl(dx p .Gl)+Gi pq Gl-Gl l Gl s ]g'" 






T ia- 




E is = 


= -vg"=E" 


w kp 


=-xG%G b sp (g ba g^+5i5" a )=-2 X G" bk G b ep =W pk 


Wis 









0<a<3, l<i,j,k,p,s<3. (dx .G j jk )=0. 



62] that the singular integral bordism groups of such an equation are trivial, i.e., 
^(JVS) ^ 0,3^ = 0. Furthermore, with respect to the notation used in diagram 
(26), one has that for the integral bordism group for smooth solutions: = 
^3w ■ Thus we can conclude that the Navier-Stokes equation is an extended 0- 
crystal PDF, but not a 0-crystal, i.e. cry(NS) ^ 0. Note that if we consider 

admissible only all the Cauchy integral manifolds X C (NS) such that all their 
integral characteristic numbers are zero, (full admissibility hypothesis), it follows 

that cry(NS) = 0. So, under this condition, (NS) becomes a 0-crystal PDE. 2S 



2 ®Let us emphasize that, similarly to the Ricci flow equation, we can identify 3-dimensional 
space-like smooth integral manifold Nq C (NS) t , such that No = Mt, via the canonical projection 
7T2 : JD 2 (W) — > M, with some smooth space-like section St : Mt C M — > Wt C W, of the 
configuration bundle it : W — > M. In fact, all the coordinates in (NS), containing time-derivatives, 
can be expressed by means of the other derivatives containing only space coordinates. (Warn for 
equation (NS)\ Even if there are not restrictions on the time derivatives of pressure, between 
space derivatives of pressure arise some constraints, as well as there are constraints between 
space derivatives of velocity components, hence sections st that identify above considered integral 
manifolds, are not arbitrary ones. For details see below.) Then we can solve the corresponding 
Cauchy problem applying Lemma 3.20, similarly to what made in the Ricci flow equation. It 
is useful to remark that in order to build envelopment solutions V = Uqgjv ^0.<j> it does not 
necessitate to handle with PDE's that admit any space-like symmetry. This of course does not 
happen for any smooth boundary value problem in (NS). With this respect, it is useful also to 
underline that the popular request to maximize entropy cannot be an enough criterion to realize a 
smooth envelopment solution. (For complementary results on variational problems constrained by 
the Navier-Stokes equation see [62].) The existence of such a smooth envelopment manifold, can 
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be proved by working on (NS) +oa . In fact, lot q, q' £ Nq C (NS) +ca , and let V q and V q i be 
two smooth solutions passing for the initial conditions q and q' respectively. We claim that their 
time-like integral curves Xo, q and X 0iq i cannot intersect for suitable short times, i.e., t £ [0, e[, 
if q ^ q'. Really if q 6 V q f| V,/ ^ 0, then T ? V g = T^V,/ = (Eoo)?. This means that in such 
a point g, Vq and V„/ must have a contact of infinity order with the Navier-Stokes equation and 

between them. We can assume that q is outside a suitable tubular neighborhood Nq x [0,e[ of 
Nq°°\ otherwise we should admit that q £ 7Vq°°'. In this last case V q and V q i should have in q 
time-like curves Xo >q and X' ^, transversal to Nq\ and with a common tangent vector C,o(q). 
Furthermore , V q and V q , should be tangent to A^ 00 ' at q: T q V q = T q V q , = (E^q D T q N^°°' > . 
Let us assume that in a suitable neighborhood JVq 00 ' x [0, e[, the manifold V q identifies two 
separated pieces, say and V^n and V 9 |2> tangent to A^ 00 ' at q and q, respectively. In other 
words VgplA^Q 00 ' X [0, e[= V,j 1 1 U ^ij 1 2 ■ (If this condition is not satisfied, then V q is necessarily 
a smooth solution of the Cauchy problem at least in a subset Y C Nq . Then we can write 

V = Vi\JV2 = {Uqgy Xo, q } LKUqgz Xo t q}, Z = A^ 00 ' \ Y, and we can continue to similarly 
discuss about the Vi-part of V.) Then, for the time-like coordinate lines ^o,g|l an d -^0,q\2t 

one has Xo l9 |i f~| -^0,g|2 = in -/Vq°°' x [0, e[. The same circumstance can be verified between 
V q i and V q , eventually by reducing e. This proves that taking e little enough, we can consider 

V = U w (oo) Xo, q C (NS),^, a fiber bundle over A^ 00 ' that is a 4-dimensional integral manifold 
representing a local solution of the Cauchy problem identified by the smooth 3-dimensional integral 
manifold No C (NS). (See Fig. 2.) On V we can recognize a natural smooth fiber bundle 
structure. In fact, since (NS) is an analytic equation and any point q 6 7Vi°°^ identifies an 
analytic solution in a suitable neighborhood U of p = ■n a a(q) £ M, it follows that if two such 
solutions s and s' are defined at q 6 Nq' they should coincide in a suitable neighborhood of p, 
denote again it by U. Then we can assume that s and s' coincide also in ft = w^(U) {~\Nq°°\ 
Therefore, we can cover 7Vq°°^ by means of a covering set {fl a }a£j , where each Q a is such that 
if q,q' S f2 Q , the corresponding analytic solutions s and s' are well defined in U a C M. In this 
way each time-like curve Xo, q is uniquely identified for q £ Q a . Then we can define smooth 
functions (transition functions) ip a p : Q a R* = R \ {0}, by r a = ipafiTp, where r a are 
nowhere vanishing sections of V — > ./Vq 00 ' on Q a C A^ 00 ' . These functions satisfy the following 
three conditions: (a) ip aa = 1; (b) Va/3 = V^t*; ( c ) ^ap^p^-ya = 1 on Q a (~| Oo (~| f2 7 (cocycle 
condition). This is enough to claim that the line bundle V —>■ JVp 00 ' has a smooth structure. Let us 
emphasize that Lemma 3.20 can be applied here also to 3-dimensional space-like smooth Cauchy 
integral manifolds N C (NS) t , that are diffcomorphic to their projections Y = n2,o(N) C Wt, 
but are not necessarily holonomic images of smooth sections of ir : W — > M. In fact, by using the 
embeddings N C (NS) C JD 2 (W) C J|(W), we can repeat above construction to build smooth 
envelopment solutions of the Cauchy problem N^°°^ C (NS) +oa C J%°(W). Really, any smooth 
3-dimcnsional space-like submanifold Y C Wt, identified by a space-like smooth integral manifold 
Z C (NS) +oa C J%°(W), such that Z = Y, via the canonical projection 7roo,o : J%°(W) -> W, 
can be locally identified by some smooth implicit functions {fi(x k , yi) = 0}i</<5, with Jacobian 
matrix of rank five, where {x , 2/ J }i<fc<3 ; i<j<5 are coordinates in the 8-dimensional affine space 
Wt, and we write Z = y(°°). Then, as a by-product we get also a cohomology criterion to 
classify envelopment solutions according to the first cohomology space H 1 (No, TL 2 ). In fact, a 
line bundle V —>■ TVq 00 ' is classified by the first Stief el- Whitney class of V, that belongs to 
H 1 (N^ cc) ,Z 2 ) ^(JVo.Za). The classifying space is ffiP 00 and the universal principal bundle 
is S°° -> RP°° = S°°/1 2 = M°°/M* S Gri(R°°), where Gri(M 00 ) denotes the Grassmannian of 
1-dimcnsional vector subspaces in R°°, and the nonzero element of Z2 acts by v 1— > —v. Since 5°° 
is contractiblc one has ir^S 00 ) S tt^RP 00 ) = 0, i > 1 and 7ri(MP°°) Z 2 - R-P°° is the Eilenberg- 
Maclane space K(X2,1) [58]. Hence [No,Gri] = H (ATq;^), / i-> where fi is the generator 
of P 1 (RP°°;Z 2 ) S Z 2 . Since one has the bijection [JV ;Gn] -> Vi(JV ), where Vi(A r ) denotes 
the set of 1-dimcnsional vector bundles over No, we get the bijection wi : Vi(A r o) — > H 1 (Nq;^) 
that is just the Stiefel- Whitney class for line bundles over Nq°°\ In conclusion an envelopment 
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In [62] it is proved that the set of full admissible Cauchy integral manifolds is not 
empty. This result gives us a general criterion to characterize global smooth solu- 
tions of the Navier-Stokes equation and completely solves the well-known problem 
on the existence of global smooth solutions of the Navier-Stokes equation. (For com- 
plementary characterizations of the Navier-Stokes equation see also [62, 67, 68, 69, 
70, 71]. There a geometric method to study stability of PDE's and their solutions, 



solution V D TV, of an admissible Cauchy integral manifolds TV C (TVS), can be identified with 
some cohomology class of H 1 (N;I>2)- In particular V is orientable iff its first Stiefel- Whitney 
cohomology class wi(V) = 0. Of course does not necessitate that all cohomology classes of 
H 1 (N; Z2) should be represented by some envelopment solution V passing for TV. However, when 
this happens we say that TV is a wholly cohomologic Cauchy manifold of (NS). This is surely 
the case when TV ~ D 3 , i.e., when the space-like, smooth, 3-dimensional, integral manifold TV 
is homotopy equivalent to the 3-dimensional disk D 3 . In fact in such a case H 1 (N;'Z2) = 0, 
hence there is an unique cohomology type of envelopment solution passing for N^°°\ (or TV), 
the orientable one. Therefore, in such a case TV is a wholly cohomologic Cauchy manifold. For 
example, when TV is identified by a smooth space-like section St : Mt C M — > Wt C W, TV is 
a wholly cohomologic Cauchy manifold. It is important to remark also that, fixed some smooth 
Cauchy problem TV C (TVS), it does not necessitate that a local smooth solution Ve should be 
unique. In fact, even if the Cartan distribution Eoo C (NS) +ao , is a 4-dimensional involutivc 
distribution, the manifold (NS), is not finite dimensional, therefore for any point q £ (NS) 
pass infinity 4-dimensional integral manifolds tangent to (Eoo)q C T q ((NS) +ao ). However, we 
can see that if there are two smooth solutions V e , V £ ' C (NS) +oc for a fixed smooth Cauchy 

problem TV C (NS), their Stiefel-Whitney classes must coincide with a same cohomology class 
of iT 1 (TV,Z2). In fact, let us denote by Noo(TV) e the integral e-normal bundle of N^°°\ i.e., 
N 00 (TV) £ = U 9eJV (oo) N 00 (TV) E , 9 , with N 00 (TV) E , 9 = ((E^q/TqN), = {[0,e[Co(g)}. Then one 
has the canonical isomorphism of line bundles over TV' 00 ': N 00 (TV) e = V e , \(o(q) i-> Xo <q (X), 
A S [0, e[, A similar isomorphism can be recognized between Noo(TV) e and V e ', taking into account 
that also Vj can be considered a line bundle: V e ' = [Jq e N^X' q . This proves that must 
necessarily be iAii(V e ) = u)\ (Noo (TV) e ) = u)\(V^). 

So we have shown that any space-like smooth 3-dimcnsional integral manifold TV C (NS) 
that is diffcomorphic to its projection on Mt C M, (for some t), via the canonical projection 
7T2 : JD' 2 (W)t — > Mt, or is diffcomorphic to its projection on Wt, via the canonical projection 
T2,o ! JD' 2 (W)t — > Wt, admits smooth solutions, V D TV. Since such integral manifolds (Cauchy 
manifolds) are not arbitrary ones, but satisfy some constraints, here we shall explicitly prove that 
such Cauchy manifolds exist. From Tab. 1 we can write (NS) in the form reported in (55) 



(55) 



From the prolongation of equation (55) (C) with respect to space coordinates x k , we get 

!( c ) 9 4 q = l p [-iq&s ~ KHi.,< - 2a^VGi - ^^ P G\ s q - .Kri,, - 
-±%Si s - x*S{ s q - ±l q T" - xlTjl - p iqg v - Pi g% 
-p(dx q dxi.f)gv - p(dxi.f)g%] 
— H ,q £ A[[x , X^ , X rs , X rsq , Pi , Piq]] 

By using equations (55)(C) and (56)(C) 9 , we can rewrite equation (55)(i?)o in the form reported 
in (57). 

(57) (E) G] k H k (x k , x a , x% , x° hk , Pl ) + H° 3 (x k ,x° ,±\ , x° hk , x s hkq , Pl , p lq ) = 0. 



(A) 


± k G j )k +±i =0 




(B) 






(B) h 




= 


(C) 


x' = i [-i s i?| - x'x^Gi^ - 
sff £ A[[x k ,x*,x* s , Pi ]] 


x s x J s p — x k S J k s — x 3 is T ls - 


(D) 


00 = ik; [-pCpx k e k - 6 is E 
= k e A[[± k ,x k ,e k ,s is ]] 


s - x k ±PW kp - x k x s p W p ks 
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Therefore, the parametric equations of a space-like integral analytic 3-dimensional manifold JV C 
(NS), diffcomorphic to Mt, has the parametric equation given in (58). 



& — (x k ) 
p — p(x k ) 

e = e\x k ) 
*l = «?*>(**) 

± =Hi(x k ) 

Pa =p (x k ) 

Pk =P,k(x k ) 

00 = K(x k ) 

e k =e k (x k ) 

(58) { 4 h =Hi h (x k ) 

±3 hk = v Mk( xk ) 

±00 = h3 o = (d±s.Hi)x' + (dii.Hi)x' h0 + (dx h a k .Hi)x' hk0 + (dp'.Hi)p i0 
= {dx s .H')H s + (d±^.W)(dx h .H") + (d±i k .H 3 )(dx k dx h .H>) 
+(dp\H')(dx i .p a ) 
Pok = (dx k .p ) 
Phk = (dx h dx k .p) 
poo = Paa(x k ) 
hk = (dx h dx k .6) 
8 0k = (dx k .K) 

OO = (dx k .K)x k + (dx p k .K)x k + (de k .K)g k0 + (d0".K)e is a 

= (dx k .K)H k + (d±l.K)x p .H k ) + (dS k .K){dx k .K) + (d$ is .K)(dXidx s .K) 

where v 1 = x J (x k ) are solutions of the continuity equation (55) (A). Therefore, equation (55) (A) 
can be considered a first order equation on the fiber bundle Et = Mt xI-> Mt, reported in (59). 



(59) 



Si C JD(Et) : {x k G 



0}. 



One can see that E\ is an involutive, formally integrablc and completely integrable PDE. In fact, 
one has 



dim(Bi) + i = 29 



= dim .Ei = 14 



dim(ji 



hence the mapping (i?i)_)_i 



dim(gi)+i = 15 



Ei is surjective. Furthermore, we get 

+ 



dim gi — 8 + 



dim^i)* 1 ) : 



dim( 9 i)( 2 > = 2 + dim(£)i)< 3 > = 



This is enough to state that g\ is an involutive symbol. Therefore, E\ is formally integrable and, 
since it is analytic it is also completely integrable. As a by-product, we get that for any analytic 
solution x k = x k (x l ) of Ei, we can write: 



(60) 



xl = Ai(x a ) - gVpi 



1 13 



where A J = A 3 (x s ) are suitable analytic functions. So equation (57) (E) can be rewritten as a 
second order equation E2 for a function p = p(x k ) as section of the trivial fiber bundle Ft = 
M[Xf-> Mf. 



(61) E2 C JD (Ft) : {g la p w +piC l {x k )-B(x k ) = Q}. 

where C and B are given analytic functions of x k . This is an involutive, formally integrable PDE, 
hence it is completely integrable, since it is analytic. In fact, 



dim(B 2 )+l = 19 



dim(B 2 ) = 12 



dim(j 2 )+l = 7 



Therefore, the map (£2)+! — > E2 is surjective. Furthermore, 



dim(p2)+i = 7 



dim((?2) — 5 



dim( S2 ) (1) = 2 + dim( ff2 ) (2) = 



hence is an involutive symbol. This concludes the proof. Thus, for any point g 6 £1, pass 
solutions of the continuity equation, and fixing a point q on the infinity prolongation (Ei)^.^, 
one identifies an analytic solution defined in a suitable neighborhood of p = Tv 00 (q) £ Mt. (Of 
course also for equation E\ we can identify smooth solutions by solving lower dimension Cauchy 
problems, by means of Lemma 3.20, i.e., by using envelopment solutions.) Similar considerations 
can be applied to the equation E2 to identify analytic and smooth solutions of the pressure 
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related to integral and quantum bordism groups of PDE 's, has been introduced, and 
applied to the Navier-Stokes equation too.) 




Figure 2. Construction of envelopment solution V for Cauchy 
problem A^ 00 ' C (NS) +oa , on the infinity prolongation of (NS) +oc 
of (NS). One has V = U qeN (°°) ^o,q, w ith X , q the unique time- 
like curve passing for q, tangent to the vector Co(<?) and contained 
in the unique analytic solution passing for q. 



4. EXTENDED CRYSTAL SINGULAR PDE's 



Singular PDE's can be considered singular submanifolds of jet-derivative spaces. 
The usual formal theory of PDE's works, instead, on smooth or analytic subman- 
ifolds. However, in many mathematical problems and physical applications, it is 
necessary to work with singular PDE's. (See, e.g., the book by Gromov [25] where 
he talks of "partial differential relations", i.e., subsets of jet-derivative spaces.) So 
it is useful to formulate a general geometric theory for such more general mathe- 
matical structures. On the other hand in order to build a formal theory of PDE's it 
is necessary to assume some regularity conditions. So a geometric theory of singular 
PDE's must in some sense weak the regularity conditions usually adopted in formal 
theory and admit existence of subsets where these regularity conditions are not 
satisfied. With this respect, and by using our formulation of geometric theory of 
PDE's and singular PDE's, we study criteria to obtain global solutions of singular 
PDE's, crossing singular points. In particular, some applications concerning sin- 
gular MHD-PDE's, encoding anisotropic incompressible nuclear plasmas dynamics, 
are given following some our recent works on this subject. The origin of singu- 
larities comes from the fact that there are two regions corresponding to different 
components PDE's having different Cartan distributions with different dimensions. 
However, by considering their natural embedding into a same PDE, we can build 
physically acceptable solutions, i.e., satisfying the second principle of the thermo- 
dynamics, and that cross the nuclear critical zone of nuclear energy production. 

functions p = p(x k ). This assures that one can identify space-like smooth Cauchy manifolds in 
(NS), that are diffcomorphic with their canonical projections on Wt, for any time t. 
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A characterization of such solutions by means of algebraic topological methods is 
given also. 

The main result of this section is Theorem 4.30 that relates singular integral bordism 
groups of singular PDE's to global solutions passing through singular points, and 
Example 4.33 that for some MHD-PDE's characterizes global solutions crossing the 
nuclear critical zone and satisfying the entropy production energy thermodynamics 
condition. 

Let us, now, resume some fundamental definitions and results on the geometry 
of PDE's in the category of commutative manifolds, emphasizing some our recent 
results on the algebraic geometry of PDE's, that allowed us to characterize singular 
PDE's. 29 

Definition 4.1 (Algebraic formulation of PDE's). Let n : W — > M be a smooth 
fiber bundle, dim = m + n, dimM = n. We denote by J k (W) the space of all fc- 
jets of submanifolds of dimension nofW and by J k (W) the fc-jet-derivative space 
of sections of it. Furthermore we denote by JD k (W) the fc-jet-derivative space for 
sections ofn. One has JD k (W) = J k (W) C J k (W). J k (W) is an open subset of 
J k (W). Let 2tfe be the sheaf of germs of differentiate functions J D k {W) —> R. It is 
a sheaf of rings, but also a sheaf of ^-modules. A subsheaf of ideals 03 a, o/2lfc that 
is also a subsheaf ofR-modules is a PDF of order k on the fiber bundle ir : W — > M . 
A regular solution of 03 A; is a section s : M — > W such that f o D k s = 0, V/ G 03a,. 
The set of integral points of 03 & (i.e., the zeros of 03 a, on JD k (W) is denoted by 
J(Q5fe). The first prolongation (03a;)+i o/Q5fe is defined as the system of order k + 1 
onW — > M , defined by the fonk,k-i and f*- 1 ', where f^ on D k+1 s(p) is defined by 
f( 1 )(D k+1 s(p)) = (dx a .(f o D k s(p))). In local coordinates (x a , y 3 , y{f) the formal 
derivative is given by f {1) (x a , y 3 , y° a ) = (dx a .f) + Y^^KkVpa^yj •/)■ The 
system is said to be involutive at an integral point q G JD k (W) if the following 
two conditions are satisfied: (i) 25fc is a regular local equation for the zeros of^k oi 
q (i.e., there are local sections F\, F t G T(U, 25fc) of*Bk on an open neighborhood 
U of q, such that the integral points of 25 fc in U are precisely the points q' for which 
Fj(q') = anddFiA- ■ -AdF t (q) ^ 0, that is F\, ■ ■ • ,F t are linearly independent atq; 
and (ii) there is a neighborhood U of q such that tt^^ ^(U) P| J^^S^+i) is a fibered 
manifold over Uf^J^k) (with projection irk+i t k)- For a system Q3fc generated 
by linearly independent Pfaffian forms 9 1 ,--- ,6 k (i.e., a Pfaffian system) this is 
equivalent to the involutiveness defined for distributions. 

Theorem 4.2 ([34]). Let *Bfc be a system defined on JD k (W), and suppose that 
Q5fe is involutive at q G J(Sfc). Then, there is a neighborhood U of q satisfying the 
following. If q G J((Q5fc)+ s ) and TTk+ s ,k{<l) is in U , then there is a regular solution 
s o/Q5fe defined on a neighborhood p = irk+ s ,-i(q) of M such that D k+S s(p) = q. 

Theorem 4.3 (Cartan-Kuraniski prolongation theorem [34, 58]). Suppose that 
there exists a sequence of integral points q^ of (Q5fe) +S ,s = 0, 1, projecting 
onto each other, TTk+s.k+s-iiq^) = 9^ s_1 \ such that: (a) (fSk)+s is a regular 
local equation for J((*B/-)_)_ S ) at qrW; and (b) there is a neighborhood of q^ 

2 ^For general informations on the geometric theory of PDE's see, e.g. ,[9, 12, 14. 22, 23, 24, 
25, 33, 34, 38, 92, 93, 94]. In particular, for singular PDE's geometry, see the book [58] and the 
recent papers [3, 75] where many boundary value problems are explicitly considered. For basic 
informations on differential topology and algebraic topology see e.g., [9, 24, 31, 39, 41, 47, 83, 88, 
87, 91, 92, 93, 94, 96, 97]. 
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in J(( ! B/ £ )+ s ) such that its projection under TTk+s ,k+s-i contains a neighborhood of 
q( s ~^ in J((5Bfe) + ( s _i)) and such that ir k + s ,k+s-i ■ -» n k + s ,k+s-i (^ (s) ) is a 
fibered manifold. Then, (23fc)_|_ s is involutive at q( s > for s large enough. 

The algebraic characterization of singular PDE's can be given by adopting the 
methods of the algebraic geometry, combined with the differential algebra. (See 
e.g., [58].) Let us go here in some details about. 

Definition 4.4. A differential ring is a ring A with a finite number n of commu- 
tating derivations di, ■ • • , d n , didj — djdi = 0, Vi, j = 1, • • • , n. A differential ideal 
is an ideal a C A which is stable by each di, i = 1, ■ • • , n. 

A differential ring (A, {dj}i<j< n ) identifies a subring (subring of constants): C = 
cst(A) = {a G A\dja = 0, Vj = 1, ••• ,n} C A. We may extend each di to a 
derivation of the full ring of fractions, Q(A), still denoted by di and such that 
di{a/r) = (rdia — adir)/r 2 , for any ^ r, a G A. 

Example 4.5. If K is a differential field with derivations d\,--- , <9 M and y k , 
k = 1 , • • • , m, are indeterminates over K , we set yfi = y k . Then the polynomial ring 
K[y]d = K[y k , k = 1 , • • • , m, \i = (1% ■ ■ ■ fi s , \fj,\ > 0], can be endowed with a struc- 
ture of differential ring by defining the formal derivations di = d{ + y k +udy k . Of 
course K[y]d is not a Noetherian ring. We write if = K[y k \k = 1, ■ • • , m; < 
\fi\ < q] and one has K(y q )d = Q(K[y q ]d). We set also K{y)d = Q(K[y]d). 

Definition 4.6. A differential subring A of a differential ring B is a subring which 
is stable under the derivations of B. Similarly we can define a differential extension 
L/K of differential fields, and such an extension is said to be finitely generated if 
one can find elements w 1 ,--- ,rf n G L such that L = K(r) ,n m ). Then the 
evaluation epimorphism is defined by K[y]d K[ 7 l]d CL, y k >-» n k . Its kernel is a 
prime differential ideal. 

Proposition 4.7. [58] Let < S >d denote the differential ideal generated by the 
subset S C A, where A is a differential ring. If A is a differential ring and a, b G A, 
then one has the following: 

(i) G< d v {ab)\\v\ < \n\ >. 

(ii) (d l a) 2r - 1 G< a r > d . 

(iii) If a is a radical differential ideal of the differential ring A and S is any subset 
of A, then a : S = {a G A\aS G a} is again a radical differential ideal of A. 30 

(iv) If a is a differential ideal of a differential ring A, then rad(a) is a differential 
ideal too. 

(v) One has the following inclusion: arad < S >dC rad < aS >d, Va G A, and for 
all subset S G A. 



30 If 

a is any ideal of A, the radical of a is the following ideal r(o) = y/a = {x S A\x n G 
for some n > 0} = rad(a). If \/a = a, then a is called radical ideal or perfect. One has also 
that r(o) is the intersection of all prime ideals p C A, containing o. In particular, the radical 
of the zero ideal < > is the nilradical, nil(A), of A, i.e.. the set of all nilpotent elements of 
A. Therefore nil(A) is the intersection of all prime ideals, (since all ideals must contain 0). One 
has also nil(A) C rad(A), where rad(A) is the ideal of A denned by intersection of all maximal 
ideals m C A. If o is a radical ideal, then A/ a is reduced, i.e., the set of its nilpotent elements is 
reduced to {0}. In particular A/nil(A) is reduced. If n : A — > A/a is the canonical projection, 
then ir~ 1 (nil(A/a)) = r(o). 
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(vi) If S and T are two subsets of a differential ring A, then 

rad < S >d .rad < T > f /C rad < ST >d~ rad < S >d Hrad < T >d . 

(vii) If S is any subset of a differential ring A, then we have: 

rad < S, ai, ■ ■ ■ ,a r >d= rad < S, a\ >d n • • • H rad < S,a r >d ■ 

Definition 4.8. A differential vector space is a vector space V over a differential 
field (K, di)i<i<„ such that are defined n homomorphisms di, i = !,••• ,n, of the 
additive group V such that: di{av) = (dia)v + a(dju), Va £ K, Vu £ V . Then we 
say that K is a differential field of definition. 

Proposition 4.9. [58] Let V be a differential vector space over a differential field 
K , with derivations di, i = 1, ■ • • , n, and let {ej}j<=/ be a basis ofV. Then the field 
of definition k of a differential subspace W C V is a differential subfield of K if it 
contains the field of definition of each d\Bi,- ■ ■ , d n ei with respect to {e^}. 

Definition 4.10. A family r\ = (rj , • • • , r] m ) of elements in a differential extension 
of the differential field K is said to be differentially algebraically independent ( or 
a family of differential indeterminates,) over K, if the kernel of the evaluation epi- 
morphism K[y]d — > K[rj\d is zero. Otherwise the family is said to be differentially 
algebraically dependent (or differentially algebraic ) over K. 

Proposition 4.11. If K/k and L/k are two given differential extensions with 
respective derivations dx and d^, there always exists a differential free composite 
field of K and L over k. 

Proof. The ring K ^ K L has a natural differential structure given by d(a ®b) = 
(dxa) <E) b + a (d^b), as 4k\k = di\ K = d. On the other hand there is a finite 
number of prime ideals pi C K(g) K L such that HiPi = an d Pi + Pj =< 1 >j 
\/i ^ j. Now we have the following lemma. 

Lemma 4.12. Ifai,--- , o r are ideals of a differential ring A such that ai + aj = A. 
Vi j, and ai PI • • • D Or is a differential ideal of A, then each ai is a differential 
ideal too. 

Therefore we can conclude that each pi is a differential ideal, hence the proposition 
is proved. □ 

Lemma 4.13. A family rj is differentially algebraic over K iff a differential poly- 
nomial PEp exists such that (dyp.P) p, where yp is the highest power of y p 
appearing in P. Sp = (dyp.P) is called the separout of P. (The initial of P is the 
coefficient of the highest power of yp appearing in P and it is denoted by Ip. More 
precisely one has P = Ip(yp) r + terms of lower degree.) 

Proposition 4.14. [58] If S is any subset of a differential ring A and r > is any 
integer, we call r-prolongation of S, the ideal 

(S) +r =< d v a\a e 5, < \u\ < r >c A. 

One has the following properties: (i) (S') + ( r+s ) = ((S)+ s )+ r . (ii) (S l )+ 00 =< S ><j- 
(iii) Let a be a differential ideal of the differential ring K [y]d- We set a q = ar\K[y q ]d, 
do = a n K [y]d, cioo = a. We call the r-prolongation of a q , the following ideal: 

(a q )+r =< d v P\P 6a„0< \v\ < r >C K[y q+r ]d- 
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One has: 

(a q ) +r C a q+r , (a g ) +oc C a, (a q ) +r n K[y q ] d = a q , Vg, r > 0. 

With algebraic sets it is better to consider radical ideals. Hence if r C K[y]d 
is a radical differential ideal, then x q is a radical ideal of K[y q ]d, for all q > 0. 
Then if E q = Z(x q ) is the algebraic set defined over K by x q = I{E q ), we call 
r -prolongation of -E g the following algebraic set: {E q ) +r = Z((x q ) +r ). In general 
one has (t g )+ r C r g + r , hence rad((r 9 )-|- r ) C x q+r . Therefore, in general one has: 

E q +r C (-^t/)+r- 

Proposition 4.15. [58] Let p C K[y]d be a prime differential ideal. Then we 
can identify each field L q = Q(K[y q ]d/p q ) with a non- differential subfield of L = 
Q(K[y]d/p) and we have: KC.LqC.---C. L^ = L. Then there are vector spaces 
R q over L q or L defined by the following linear system: 

(dy%.P T )(r])v*=0, {l<r<t, l<fc<m, \fi\ = q} , 

where i] is a generic solution of p and Pi,-- - , Pt are generating p q . Such result 
does not depend on the generating polynomials. We can also define the vector space 
g q (symbolj over L q or L, by means of the linear system: 

{dy£.P T )(ri)v* = 0, {1 < r < t, l<k< m, < H < q} . 

For the prolongations (g q )+ r one has, in general, g q + r C (g q )+ r , Vq> f > 0. 

Definition 4.16. We say that R q or g q is generic over E q , if one can find a 
certain number of maximum rank determinants D a that cannot be all zero at a 
generic solution of p. 

Proposition 4.17. R q or g q is generic if we may find polynomials A a , B T £ K[y q ]d 
such that: 

Y J A a D a + Y J B T P T = 1. 

a r 

Furthermore, R q or g q are projective modules over the ring K[y q ]d/p q C K[y]d/p. 

Proof. It follows directly from the Hilbert theorem of zeros. (See [58].) □ 

Theorem 4.18 (Primality criterion [58]). Let p q C K[y q ]d and p 9 +i C K[y q+ i] d 
be prime ideals such that p g +i = {p q )+i and p q +i f) K[y q ]d = p q . Lf the symbol g q of 
the variety R q defined by p q is 2-acyclic and its first prolongation g q +\ is generic 
over E q , then p = (p q )+oo is a prime differential ideal with p n I"([y q -\- r ]d = (p q )+ r; 
for all r > 0. 

Let x q C K[y q ]d and r 9 +i C K[y q+ i]d be radical ideals such that t 9 +i = (r)+i and 
t g +i C\K[y q ]d = Xq. If the symbol g q of the algebraic set E q defined by x q is 2-acyclic 
and its first prolongation g q +i is generic over E q , then x = (r 9 ) +00 is a radical 
differential ideal with x D K[y q+r ]d = (x q )+ r , for all r > 0. 

Theorem 4.19 (Differential basis). If x is a differential ideal of K[y]d, then x = 
rac?((r g )-)- 00 ) for q sufficiently large. 

Proof. In fact one has the following lemma. 

Lemma 4.20. If p is a prime ideal of K[y]d, then for q sufficiently large, there is 
a polynomial P £ K[y q ]d such that P ^ p q and Pp q + r C rad((p q )+ r ) C p q + r , for 
all r > 0. 
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After above lemma the proof follows directly. □ 

Every radical differential ideal of K[y] can be expressed in a unique way as the non- 
redundant intersection of a finite number of prime differential ideals. The smallest 
field of definition k of a prime differential ideal p C K[y] is a finitely generated 
differential extension of Q. 

Example 4.21. With n — 2, m = 1, q = \. Let us consider the differential 
polynomial P = y\y\ — y\y\ ~ 1- We obtain for the symbol g\: y\v\ + y\v\ — 
y\v\ — y\v\ = 0. Setting v f = yfw\ we obtain (y\y\ — y\y%){w\ + w\) = and thus 
w 2 + w l = on Ei. Hence gi is generic. One can also set Pi = y\, P2 — y\ and 
we get the relation: y\Pi — y\P2 — P = 1. A similar result should hold for E\. gi 
is involutive and the differential ideal generated by P in Q < y l ,y 2 > is therefore 
a prime ideal. 

Definition 4.22. A differentially algebraic extension L over of a differential field 
K is a differential extension over K where every element of L is differentially 
algebraic over K . 

The differential transcendence degree of a differential extension L/K is the number 
of elements of a maximal subset S of elements of L that are differentially transcen- 
dental over K and such that L becomes differentially algebraic over K(S). We shall 
denote such number by trdd{L/K). 

Theorem 4.23 ([58]). One has the following formula: 

. ^ ( r _|_ {\\ 

dim(p 9+r ) = dim(p 9 _i) + ^ — - a l q , Vr > 0, 

l<i<n 

where a* is the character of the corresponding system of PDE's. The character a l q 
of a q-order PDE E q C JD q (W), n : W — )• M, dimM = n, with symbol g q , is the 

integer a\ = dim(5^ _1) ) p - dim(.g^ ) ) P; p € E q , where {g q l) ) p = {C G i,9q)p\Q{vi) 
■ - ■ Ci v i) = 0}, where (vi, ■ ■ ■ ,V n ) is the natural basis in T„ k ^M . 
The character a™ and the smallest non-zero character only depend on the differential 
extension L/K and not on the generators. In particular, one has: trdd{L / K) = a q . 
If C is differentially algebraic over K{rf}d and rj is differentially algebraic over K, 
then Q is differentially algebraic over K . 

If L/K is a differential extension and £,77 G L are both differentially algebraic over 
K, then £ + r\, £77, £,/rj, (rj ^ §), and di^ are differentially algebraic over K . 

Theorem 4.24. Let [A, {dj}i<j<n) be a differential ring. The set D(A) of dif- 
ferential operators over (A, {dj}i<j< n ) is a non-commutative filtered ring and a 
filtered bimodule over A. 

Proof. If y is a differential indeterminate over A, we may introduce the formal 
derivatives di, • ■ ■ ,d n which are such that didj — djdi = 0, Vi, j = 1, • • • , n, and 
are defined by: d.^ay) = {dia)y + a(diy). We shall write d t y = y it d l y l _ l = y^+u, 
where fi is the multi-index fj, = (fj,i, • • • , /j„) with length = fii + ■ ■ ■ + fj. n . If 
y = (y 1 , • • • ,y m ), we set d^ = {diY 1 ■ ■ ■ (d n )^ n and d^y k = y^. Any differential 
operator of order q over A can be written in the form P = J2o<fj,<q aAl ^M' flM G A. 
Set ord(P) = q. Then, we can write D(A) = A[d\,--- ,d n ] = A[d] the ring of 
partial differential operators over A with derivatives di, ■ • ■ ,d n . The addition rule 
is clear. The multiplication rule comes from the Leibniz formula: 
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Table 2. Examples of singular PDE's defined by differential polynomials. 



Name 


Singular PDE 


PDE with node and triple point 
El C JD(E) 


Pi = 04)* + ( U g)4_( u i)2 =0 
P2 = (u 2 ) 6 + K) 6 - u2uj = 


PDE with cusp and tacnode 
E l C JD(E) 


yi = K) 4 + (n2)4-(^)3 + («2)2=0 

?2 EE (« 2 ) 4 + «) 4 - («2) 2 K) - (^)K) 2 = 


PDE with conical double point, 
double line and pinch point 
Ei C JD(F) 


ri EE - («i)K)2 = 

r 2 EE(« 2 ) 2 -(n2)2-(«l)2=0 
r 3 EE(« 3 ) 3 + ( U 3)3 + (u 2 )(u 3 )=0 



vr : B EE R 4 -> R 2 , (x, j/.n 1 ,^ 2 ) i-s> vr : F EE R 5 -> R 2 , (a;, y, u 1 , u 2 , u 3 ) M- (i, jf). 

a =< pi, pa >C A, fa ee< <ji, (J2 >C A, c ee< ri, r2, r-3 >C B. 



Here we have put /i! = /ii! • • ■ fi n \. With these rules D (A) becomes a non-commutative 
ring and a bimodule over A. In fact, the previous formula defines the right action 
of A on D(A). The left action of A on D (A) is simply the multiplication on the 
left by A, that is aP = a(Eo<^<</ flMc W = Eo</i<9 aa ''^- Now, the filtration of 
D(A) is naturally induced by filtration of spaces of differential operators. More 
precisely D q (A) = {?£ D(A)\ord(P) < <?}, where ord(P) = sup{\p,\\a^ ^ 0}. We 
set D- X (A) = and D (A) = A. Then, D q (A) C D q+1 {A), D(A) = \J q > Q D q (A) 
and Dg(A)D p (A) C £> p+g (A). □ 

Theorem 4.25 (Algebraic criterion for formal integrability [58]). Let Z q = Z(p q ) 
be the variety defined by means of ideal p q C K[y q ]d such that the following condi- 
tions are verified: 

(i) (p q )+i = pq+i C K[y q+ i]d is also a prime ideal. 

(ii) pg+i n K[y g ] d = p q . 

(iii) g q +i is generic over E q . 

(iv) g q is 2-acyclic. 

Then (p 9 )+oo = P C K[y\d is a prime differential ideal, where p is the differential 
ideal generated by a finite number of differential polynomials P\, ■ ■ ■ ,Pt, defining 
E q , and E q is formally integrable. If one of these conditions is not satisfied we get 
that p is not a prime ideal, hence we have a factorization of p. In other words the 
PDE is not formally integrable. 

Proof. For a detailed proof see [58] . □ 



Example 4.26 (Some singular PDE's). In Tab. 2 we report some singular PDE's 
having some algebraic singularities. These are singular PDE's of first order defined 
on JD(E) = R 8 for the first two and on JD(F) = R 11 for the third. To the ideals 
a =< Px,P2 >C A, b =< qi,Q2 >C A and c =< ri,r 2 ,r 3 >C B, where A = 
Kfu 1 , it 2 , u x , Uy, u 2 ,, Uy] and B = M^u 1 , u 2 , u 3 , u x , u y , u 2 x , w 2 , u%, Uy], one associates 
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the corresponding algebraic sets 

{ E l = { q e M. s \f(q) =0,V/eo}c R 8 

Ei = {qe M 8 \f{q) = 0, V/ G fa} G R 8 
Ei = {q G K u |/(g) = 0,V/ 6c}c R 11 . 



(62) 



These algebraic sets are in bijective correspondence with the corresponding radicals: 
r{a) = {g G A\g{p) = 0,Vp G E x } D a, r(b) = {g G A\g{p) = 0,Vp € Ex} D fa,. 
r ( c ) = {9 6 B\g(p) = 0,Vp G -Ei} D c. (This follows from the Hilbert theorem of 
zeros [58] .) 

Let us consider in some details the singular PDE E\ C JD(F), in order to see 
existence of global algebraic singular solutions and characterize their stability. We 
have the following representation: E\ = (J A\ (J A 2 [J A 3 , where 



(63) Aj 



I) 



2) u\ = «(j)V(« 2 ) 2 - ( U D 2 



3) ul 



ju 3 ) 3 



«) 



3\2 



> C Ex G JD{F) 



(64) A 3 = {q G JD(F) S R 11 ^ 1 = u 2 = u 3 = u 2 x = u\ = u 2 y = u 3 = 0} . 

In (63) we put s(l) = 1 and s(2) = — 1. A3 = R 4 is the set of singular points of 
Ei . Instead A\ and A2 are formally integrable and completely integrable PDE 's. In 
fact, one has the exact commutative diagrams (65) and (??), with j = 1,2; 

(65) 



A 3 



A 3 



M 



JD(F)^F -0 JD(F)^M 



This can be seen rewriting equations in (63) in the following equivalent way: 
(66) 



n v l - {u ).. 

1) 



Ai={ 2) ul = s{j)J{u*y-[^- + {ul) 



3) u»=-[^ + («J) 2 ] 



(ulf < (u 2 ) 2 

M 3 ^0 



1 C-Ei C JD{F). 
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The first prolongation (Aj)+i, j = 1,2, is given by the following equations: 
(67) 
{ 1) n = 

2) r 2 = 

3) r 3 = 



4) u\ x = pu^K) 2 " (^^^/(u 2 ,) 4 

5) < = [2«XK) 2 - MW^l/K) 4 



6) "xy 
7) 



where 



s(i)f/^ 2 ) 2 -[^ + K) 2 ] 2 



(« 2 ) 2 < (^ 2 ) 2 



« 2 ^o 



9) <, 



K) 2 



I 



+ 2ulu 3 yx ] 



c J£ 2 (F) 



,a..a or(« 3 ) 3 , /--3>2 ir 3 C« 3 ) 3 "l»g-(" 3 ) 3 "g, 

K) 2 



a = 2^-2[^ + (^) 2 



h = 2u>ul - 2[*j£ + K) 2 ][ 3( " 3)2( t;^ )3 "- + 2^ 



Tften 6j/ using (67) (9) to substitute u yx in (67) (1), and by using the two different 
expressions of u\ y in (67)(5) and (67)(6) to obtain an explicit expression of u yy , 
we get the following equations for the first prolongation of (Aj) + i: 



(69) 



1) ri=0 

2) r 2 = 

3) r 3 = 



4) u\ x = \2v>u\{u%f - (uTa^l/K) 4 

5) ^^Wf/^ 3 ) 2 -^ + K) 2 ] 2 



with 
(70) 



• 6 ) «™ = ^)l/^(« 2 ) 2 -[^r + K) 2 ] 2 

r3(« 3 ) 2 ^^-(» 3 ) 3 ^ , . ,3, i 



7) «L 

8) « 2 



r 3(« 3 ) 2 K) 2 -(« 3 ) 3 « 3 



fem + Au y u yy\ 



K) 2 < (« 2 ) 2 



[ ^ VJL +^ yUyy ] 



= (Aj) +1 c JD 2 (F) 



EXTENDED CRYSTAL PDE'S 



41 



Therefore also (Aj)+i are analytic submanifolds of JD 2 (F), for j = 1,2. Further- 
more, since (dim(Aj)+i = 11) = (dirr^A, = 8) + (dim(<?i)+i = 3), we see that 
the canonical projections 7T2,i : (Aj)+i — > Aj, j = 1,2, are affine subbundles of 
JD 2 (F) — > JD(F), with associated vector bundles (<7i)+i — > Aj, j = 1,2. Fi- 
nally the symbol dim(<7i) ggj 4 = 3, j = 1,2, and dim(dx\ (gi)q<=Aj ) = 0. Therefore, 
(9i)q£Aj o,re involutive. This is enough to conclude that Aj are formally integrable, 
and since they are also analytic they are completely integrable too. Furthermore, 
taking into account that dim A, = 8>2x2 + l = 5, we can apply Theorem 2.15 in 
[60], (here reported in Section 3), to calculate the weak and singular integral bor- 
dism groups of Aj. One has n i 3 w = fl 1 3 s ^ 0, j = 1, 2. Therefore, Aj are extended 
0-crystal PDE's. So for any two admissible closed 1-dimensional smooth integral 
manifolds No,N% C Aj, there exists a (singular) 1-dimensional integral manifold, 

solution V C Aj , such that dV = Nq (J N\. 31 Such a solution is smooth iff all the 

integral characteristic numbers of Nq are equal to ones of Ni . 

The Cartan distribution on E\ is given by the following vector fields 

(71) C - X x {8x + u k x du k ) + X y {dy + u k y du k ) + Y k x du k + Y k v du k y 

such that the following equations are satisfied: 

{ Yf(u 2 y f + y a w 2uiuj - 2u\ulx* + u\xy) = 



(72) 



y?2u£ + Y 2 x 2u 2 x - 2u 2 {u 2 x X x + u 2 y Xy) = 
. Y 2 X u 3 y + Y 3 y (u 2 x + 3(u 3 y ) 2 ) + 3{u 3 ) 2 (u 3 x X x + u 3 y Xy) = 0. 



Therefore dim(Ex) g eA.j = 5, j = 1,2. For example on A\, one has the following 
expression of the Cartan vector field: 



(73) 

' C = X*[dx + u k du k + ^dul + ^<du\ - jSS^dul] 
< +XV[dy + u k du k + 2 ^dul + ^du\ - [iSSf^] 

+Yf [du 2 x -Otdul- [ul+ % iyi] du y ] + Y 3 x du 3 x + Yi [- 2 -<dul + du 2 y ] . 

Instead, since equations (72) are identities for q € A 3 , we get that (Ei) 9e ^ 3 = 
~Ei(F) q £A a , i.e., in the singular points the Cartan distribution of Ei coincides with 
the Cartan distribution of JD(F) that is just given by vector fields given in (71) 
for arbitrary functions X x , X v , Y k x , Yi : JD(F) — > M, k = 1,2,3. Thus we can 

prolong any solution V C Aj, j = 1,2, dV = Nq[JNi to a solution Z, such that 
dZ = N% U{qo} = Ni, where go S A3 also. In other words, V' = V[J Ni Z is 
an algebraic singular solution of E\. (See Fig. 3.) In fact, we can always find 
a solution V C J^F) of the trivial equation J\(F) C J^F), such that dV = 



J Let us emphasize that by using Lemma 3.20 we can identify admissible smooth 1-dimensional 
integral manifolds in Aj, j = 1, 2. In fact, Aj are formally integrable and completely integrable 
PDE's. So we can use Lemma 3.20(1), but also Lemma 3.20(3), since irifi^Aj) = F, and Aj — > F 
are affine subbundles of JD(F) — > F, with associated vector bundle the symbol g\. 
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Figure 3. Algebraic singular solution V' — ^Ujvi 
JD(F), passing through a singular point q G A 3 . 



Z C E 1 C 



Ni[J{qo} = N±, and such that there exists a disk D\ C V, centered on qo, with 
radius e and boundary dD\ = N e C Aj. (Let us emphasize that dimF = 5, hence 
we can embed in F any 2-dimensional smooth compact manifold. See, e.g., [31].) 

Let V be the submanifold of V such that dV = N e {JNi. Then, since Aj is a 
strong retract of J^F), we can deform V, obtaining a solution V" C Aj, such 
that dV" — N C [J N\. By taking the limit e — > 0, we can see that the solution V" 
identifies a solution V" C Aj, such that V"{J {qo} = V C E\ is just an algebraic 
singular solution of the singular equation E\ . 

Note that the symbol (g\) q ^A 3 = G\ = TM®F, i.e., (gi) qe A 3 coincides with the 
symbol of the trivial PDE JD(F) C JD(F). In fact the components of the symbol 
on Ei, must satisfy the following equations: 



(74) 



Y?2u\ + Yi2u 2 x = 0; 
I Y 2 *ul + Yg[ul + 3(ulf} = 0. 



These equations are identities on q G A3. As a by-product, we get that above 
considered algebraic singular solutions are, in general, unstable in finite times. 

Definition 4.27. We define extended crystal singular PDE, a singular PDE E^ C 
J^(W) that splits in irreducible components Ai, i.e., Ek = Ui^i? where each Ai 
is an extended crystal PDE. Similarly we define extended 0-crystal singular PDE, 
( resp. 0-crystal singular PDE ), an extended crystal singular PDE where each com- 
ponent Ai is an extended 0-crystal PDE, (resp. 0-crystal PDE). 

Definition 4.28 (Algebraic singular solutions of singular PDE's). Let Ek C J*(W) 
be a singular PDE, that splits in irreducible components Ai, i.e., Ek = 1J 4 Ai. Then, 
we say that Ek admits an algebraic singular solution V C Ek, if V p| A r = V r is 
a solution (in the usual sense) in A r for at least two different components A r , 
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say Ai, Aj , i ^ j, and such that one of following conditions are satisfied: (a) 
(ij)Ek = Aif]Aj ^ 0; (b) ^'E k = Ai\JAj is a connected set, and Uj)E k = 0- 
Then we say that the algebraic singular solution V is in the case (a), weak, singular 
or smooth, if it is so with respect to the equation (;j\E k . In the case (b), we can 
distinguish the following situations: (weak solution): There is a discontinuity in V , 
passing from Vi to V}; (singular solution): there is not discontinuity in V, but the 
corresponding tangent spaces TVi and TVj do not belong to a same n-dimensional 
Cartan sub- distribution of J*(W), or alternatively TVi o,nd TVj belong to a same 
n-dimensional Cartan sub- distribution of J*(W), but the kernel of the canonical 
projection (717^0)* : TJ^(W) — > TW, restricted to V is larger than zero; (smooth 
solution): there is not discontinuity in V and the tangent spaces TVi o-nd TVj 
belong to a same n-dimensional Cartan sub- distribution of J„(W) that projects 
diffeomorphically on W via the canonical projection (71^0)* : TJ^(W) — > TW . 
Then we say that a solution passing through a critical zone bifurcate. 32 

Definition 4.29 (integral bordism for singular PDE's). Let Ni,N 2 C E k C J*{W) 
be two (n — 1) -dimensional admissible closed integral manifolds. We say that N\ 
algebraic integral bords with N 2 , if N\ and N 2 belong to two different irreducible 
components, say N\ C Ai, N 2 C Aj, i ^ j, such that there exists an algebraic 

singular solution V C Ek with dV = Ni \J N 2 . 

In the integral bordism group fi^-i ( res P- ^n-i s> res P- ^n-i w) °f a singular PDE 
Ek C J*(W), we call algebraic class a class [N] G ^n-ir ( res P- [N] £ ^n-i. s > resp. 
[N] G ^n-i>)> N C Aj, such that there exists a closed (n — 1)- dimensional 
admissible integral manifolds X C Ai C Ek, algebraic integral hording with N, 
i.e., there exists a smooth (resp. singular, resp. weak) algebraic singular solution 

V C E k , with dV = N\JX. 

Theorem 4.30 (Singular integral bordism group of singular PDE). Let Ek = 

[J i Ai C J^iW) be a singular PDE. Then under suitable conditions, algebraic sin- 
gular solutions integrability conditions, we can find (smooth) algebraic singular 
solutions bording assigned admissible closed smooth (n — 1)- dimensional integral 
manifolds No and Ni contained in some component Ai and Aj, i ^ j. 

Proof. In fact, we have the following lemmas. 

Lemma 4.31. Let E k = \J l Ai C J*{W) be a singular PDE with {ij)E k = Ai f| Aj ^ 
0. Let us assume that A { C J*(W), Aj C J*(W) and (ij)E k C J%(W) be 
formally integrable and completely integrable PDE's such that dimA^ > 2n + 1, 
dim A, > 2n + 1, dim ^Ek > 2n+ 1. Then, one has the following isomorphisms: 

\ to ) "n— 1,M> — iL n-l,w — iL n-l,w 

So we can find a weak algebraic singular solution V C Eh such that dV = No [J N\, 
with N a cA l ,N 1 c A 3 , iff Ni G [No]. 

Furthermore, if g k (Ai) ^ 0, g k+1 (Ai) ^ 0, gk(Aj) ^ 0, g k +i(Aj) ^ 0, g k ((ij)E k ) ^ 
0, g k +i((ij)E k ) ^ 0, then one has also the following isomorphisms: 

(76) n A \ ^n A > snw j f fc . 

V ^ ) n—l,s n—l,s n— 1,5 

32 Note that the bifurcation does not necessarily imply that the tangent planes in the points 
of Vij C V to the components Vi and Vj, should be different. 
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So we can find a singular algebraic singular solution V C Ek such that dV = 
NoUNi, N c A, Ni c Aj, iff ^ e [N ]. 

Proof. In fact, under the previous hypotheses one has that we can apply Theorem 
3.16 to each component A;, Aj and (ij)E k to state that all their weak integral 
bordism groups of dimension (n — 1) are isomorphic to ® r+s= „_ 1 H r (W; Z 2 )<g>z 2 ^s- 
Furthermore, under the above hypotheses on nontriviality of symbols, we can ap- 
ply Theorem 2.1 in [55]. So we can state that weak integral bordism groups are 
isomorphic to the corresponding singular ones. □ 

Lemma 4.32. Let E k = |J- A be a 0-crystal singular PDE. Let ^ E k = Ai\JAj 
be connected, and (ij)E k = Ai C\ Aj ^ 0. Then O n _^J = 0. 33 

Proof. In fact, let Y C (ij)Ek be an admissible closed (n — l)-dimcnsional closed 
integral manifold, then there exists a smooth solution Vj C Ai such that dVi = 

N \JY and a solution Vj C Aj such that dVj =Y\JNi. Then, V = V t \J Y Vj is 
an algebraic singular solution of E k . This solution is singular in general. □ 

After above lemmas the proof of the theorem can be considered done besides the 
algebraic singular solutions intcgrability conditions. □ 

Example 4.33 (Extended crystal singular MHD-PDE's). In the recent paper [68] 
we introduced a new PDE of the magnetohydrodynamics encoding the dynamic 
of anisotropic incompressible plasmas able to describe nuclear energy production. 
This equation is denoted (MHD) C JD 2 (W) and it is reported in Tab. 3. The 
fiber bundle there considered is 7T : W — > M, over the Galilean space-time M, 
with W e M x I x S 3 x R 4 ; where I is an affine 3- dimensional space (time-like 
flow velocities space) and S is a 3- dimensional Euclidean vector space. A section 
s = (v,p,6,E,H,I,p,h) represents flow-velocity, isobaric pressure, temperature, 
electric vector field, magnetic vector field, electric current density, electric charge 
density, nuclear energy production density. In that paper it is also proved that equa- 
tion (MHD) is formally integrable, completely integrable and an extended 0-crystal. 
Now, from Theorem 3.16 we can see that for any two space-like admissible Cauchy 
hypersurfaces Nq C (MHD) to , N\ C (MHD) ti , t\ ^ t2, there exists a (singular) 

solution V C (MHD), passing through Nq and N\. The admissibility requires that 
Nq and Ni are smooth 3- dimensional regular manifolds with respect to the embed- 
ding (MHD) C Jf(W), locally satisfying the Cauchy problem and with orientable 
boundaries Xq = 8Nq, X\ = dN\, bording by means of a suitable 3- dimensional 
time-like integral manifold P. 

Then the solution V has boundary dV = Nq[J Xo P[J Xi Ni. (For details see 
Refs.[69, 71] There can be also found the explicit expressions of the differential 

polynomials defining (MHD).) In [71] we proved that we can identify a sub- 
equation (AIHD) C (M HD) such that in some neighborhood of its points, there 
exist entropy-regular-solutions passing from such initial conditions. More precisely, 



'But, in general, it is Q n _i k ^ 0. 
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Table 3. Nuclear energy producing magnetohydrodynamics equa- 
tion {F^ = 0; 1 < s < 7} : (MHD) c JD 2 (W). 



Maxwell 



F' 1 ) = F^ 1 ' = B k / k = (no magnetic monopoles) 

D k /k ~ P~47r = (Gauss's law of electrostatic) 



= F (2) 
p(3) = F< 3 )' 

F( 4 ) = F< 4 )' 



e ijs E j/s + \{dt- B i ) = (Faraday's law) 



4?r j 



(Ampere's law) 



D k 

r t k 



kt Ei, (e =dielectric permeability tensor) 

- - -eo9 ife + e( v r/s + v s/r)9 rl 9 sk 
B k = [i kl Hi, ((j, =magnetic permeability tensor) 



i k 



Navier-Stokcs 



v n /k = (continuity equation) 



i?(5) = p(5) = 

(first-prolonged continuity equation) 



F(r) 



' jk — F(t, oc iy) = ®> (motion equation) 



— - »(0/i)/k9 ik + S k - [2 X e lk + j^{B i B k 



P Cv It T "It 
+PE : > gij — ph = 0, (energy equation) 
<j fc = —u8/ i g lk , (heat flow) 



■F i F fe )]i; i 



/fc 



r/s+v 3/r )g"g° k +M 



tk 



ki 



lk (full stress tensor) 
-^{B^-Bi +E l E3) (deviatory stress) 



x{v r / 3 +v a / r )g rl g : 
(magnetic stress tensor) ( B )jtf« = \B k B k g 13 + B l B 3 ) 

(electric stress tensor) ( E ~IM^ = -^(-\E k E k g ij +E i E^>) 
(e.m. stress tensor) Af« = ( B ~> APJ + (^M™. 



(body force): F, 



(body) 



-p{dx k .f)g ki + pE* + e i i k I J B k 



h body energy source density. 



we can represent 9 2 1Z as a polynomial differential of first order on JD(W), belong- 
ing to A[x a ,x%, H p , E p ,9,9j,h] ; where A = [[x , x 2 , x]]. This can be seen taking 
into account that 

(77) 9 2 K(s) = 9 -^-{B l Bi + E'E 1 )^ + —<"■'<,, + h + -(gradfl) 2 . 

and by calculating the corresponding explicit differential polynomial expression for 

9 2 K : JD(W) -> R. (See Tab. I) (MHD) C (MHD) is an extended 0-crystal 
singular PDE, in the sense of Definition 4-27. 

Let us resume the proof given in [71], since this is necessary to understand the 
further developments reported below. Let us define 

(78) (MHD) = jg g (MHD)\9(q) > 0, h(q) > 0, K(q) > j . 

Then, (MHD) is a connected, simply connected bounded domain in (MHD). In 
fact, we can split (MHD) in the following way 

(79) (MHD) = (+,+)?! |J (+i+) Y 2 |J Y 3 |J (+ ^ +) (MHD) , 
where 

(80) (+:+ . +) (MHD) = { g e (MHD)\h(q) > 0,9(q) > 0,1l(q) > o} 
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Table 4. Polynomial differential expression for 6 2 1Z. 



e 2 n = e [x a b x d n(i) bd c + xix d wwn{2) b a % q + x%x<>K(3) b c + xix-Hvmn{i) b acpq 

+x a x b K(5) ab + x a x b Tl{&) abpq + xfH°H r 1l{7)mr + x™E*Em(8y ml:j 

+x t HiH r n{V)t q r+x t E i E31l(10)ti 3 + x r u xlx d HPHiTZ(ll)^ pq + ±r u xl± c HP HiK{12)y b cpq 
+x r u x a x b HPHiTl(n)]i abpq + x m x a x b HPH"TZ(14) vjabpq + fej + 9^11(15)^ 

ni) b a d c = ^(9 bd 9ac + S d J b ) 

KMaipc = -^(SiS b g pc + 5 b 5^g ca + g db g qa g P c) 
^(3)L = f G b ac 

^■(^)acpq = ~^^(^q9pkGac + ^q9k a G pc + S p g a jG qc + 2g pa G qc + 2g a qG pc + gqtgajG^G^) 

K(5) ab = ?fG\ a Gl b 

K(6) abpq = -tm(g kjGqa G k b + g qt Gi a G k b + 2g pk G j qa G k b ) 

2 

K(8); ;j = sh(S s jg i m + S!g jm ) 

n{9) tqr = ^ p g m rG™ 

niO) U] = s^(g im Gni + g jm G%) 

n^)t h a d c V q = £^(S%S b 9 gic + SXg ca + g db gqag lc W r S; + g™g rp ) 

K(12)? b opq = £-[±(S b g ik G k c + S\g ka G k c + 8 b g aj G qc + 2g ia G b c + 2g aq G b c + gqtg a jG kc G k b ) 

+(6%S b g ir + S b 5^g ra + g ab gq a g, r )G pc ] 
K(13)? abpq = -§^p-[\{G qa G k b g k j + G ka G k b g qt + 2G{ a G k b g ik ){& i v 6 p +g iu g rp ) 

+^g lk G k b + S»g kr G k b + ^g r] G ] qb + 2g lr G^ qb + 2g rq G k b + g q tg rj Gl b G k u )G pa ] 
^(14)„ai, P9 = T^(G\ a G k b g ki + G\ a G k b gqt + 2G qa G jb g ik +)G pw 
TZ(15y k = -g kj 



is an open submanifold of (MHD), and 

\ (+,+)?! ={ge(M^D)|%) = 0,%)>0,^(g)>0} 

(81) I ^ 

[ {+,+)% (MHD)\1Z(q) = 0, 6(g) > 0, h(q) > o} 

are codimension 1 submanifolds of (MHD) . Furthermore, 

(82) % = [qe {+) (MHD)\h(q)^Q 1 'R(q)=Q) 
is a codimension 2 submanifold of 

{+) (MHD) = \qe (MHD)\0(q) > o} . 

Let us study the integrability properties of such submanifolds of (MHD). First note 

that since ( + . +j+ )(A/iJZ?) is an open submanifold in (MHD) and in JD 2 (( +i+ ) W), 
where 

(+,+)W = M x I x R x R+ x S 3 x R x R+, 

it follows that if q 6 ^ +i++ ^(MHD) there exist solutions belonging to ( + ^ +}+ )(MHD), 
i.e., entropy-regular-solutions passing through q, just follows from the fact that 
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(+,+,+){MHD) is an open submanifold of (MHD). So ( +>+!+ \(MHD) is com- 
pletely integrable. Furthermore, since 



(+,+,+)] 



(MHD) +r = JD r ( {+ ^ +) (MHD) ) f| JD^{ [+ ^W) 



(83) 

C JD r ((MHD))f]JD 2+r (W) = (MHD)_ 
it follows that ^ +:++ ^(MHD) is also formally integrable. Let us consider the 
integrability properties of (+.+)Yi C d(MHD). If q G (+ t +)Yi it follows that 
q G ( + }{MHD). Since this last equation is formally integrable and completely 
integrable, and ( +i+ )Yi is an open submanifold of (MHD), it follows that also 
r+ t+ \Yl is formally integrable and completely integrable. 
Let us, now, study the integrability properties of 

(84) (+>+) y 2 C JD 2 ( ( ^ +) W) :{ll = 0, = o} . 

One can see that this equation is not formally integrable, but becomes so if we add 
the first prolongation of1Z = 0. So we can prove that the following equation 

(85) (+i+) y a g jd 2 ( {+ , +) w) -.[n = Q,n a = 0, F^ = o} 

is formally integrable and completely integrable. Similar considerations can be made 

on the part Y$, and we can identify, the corresponding formally integrable PDF Y3. 
We skip on the details. By conclusion we get that 

(86) (MHD) ee (J (+ , +) y 2 |J % |J (+i+t+) (Mgg) c (MHD) 

is the formally integrable and completely integrable constraint in (MHD), where 
for any initial condition, passes an entropy-regular-solution for (MHD). Then we 
can apply Theorem 3.16 to conclude that (MHD) is an extended crystal singular 
PDF. Furthermore, taking into account that it results 

we get that (MHD) is also an extended O-crystal singular PDE. This assures that 

for any two space-like admissible Cauchy integral manifolds N\ G (MHD) , N2 G 

fl 

(MHD) , t\ 7^ t2, such that if both Ni, i = 1,2, belong to the same component, 

in the split given in (86), there exist (singular) entropy-regular-solutions V such 
that dV = Ni[J Xi P[J X2 N2, and such that their boundaries X\ ee dN\, X2 = 
dN2, should be orientable and propagating with an admissible 3- dimensional integral 
manifold P. (For details see [69]. ) The stability properties of (AIHD) and its 
solutions, can be studied by utilizing our recent geometric theory on the stability 



34 Let us emphasize that admissible Cauchy integral manifolds can be found in each component 
of (MHD), thanks to Lemma 3.20. More precisely, the proceeding followed for the Navier-Stokes 

equation in Example 3.27, to solve Cauchy problems there, can be applied also to (MHD). In 
fact, envelopment solutions can be built also for all the components of this last equation, since 
they are formally integrable and completely integrable PDE's. 
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of PDE's [67, 68, 69, 70, 71]. More precisely, (MHD) is a functionally stable 
singular PDE, in the sense that it slpits in components that are functionally stable 
PDE's. Furthermore, smooth entropy-regular-solutions, i.e., smooth solutions of 

(MHD), do not necessitate to be stable. However, all they can be stabilized and 
the stable extended crystal singular PDE of (MHD), i.e., a singular PDE slpits in 
components that are stable extended crystal PDE's. This last is just ^(MHD) = 
(MHD) . There, all smooth entropy-regular-solutions, belonging to only one of 
the components, in the split representation (86) are stable at finite times. 
Let us, now, find global solutions of (MHD), crossing the singular sets correspond- 
ing to states without nuclear energy production, to ones where h > 0, i.e., passing 

from ( +j+ )yi to ( + + + -)(MHD). Let V C (MHD) be a time-like solution such that 
the following conditions are satisfied: 

(i) V is a regular entropy solution; 

(ii) Vr\(+,+)Yi= {+) V^<D; 

(iii) V f] (+,+,+) (MHD) ee (+t+) V ± 0. 

(iv) dV = N {jP\jN l , N c (+) V, N\ C (+)+) V, d (+) V = N \J (+) P\JY. 
Then we can give the following split surgery representation of V : 



(87) V = (+) V{JZ 

Y 

with Z = Y (J (+,+) V. We call a solution V of (MHD), such that holds the surgery 
property given in (87) just a crossing nuclear critical zone solution. Then, for any 

two admissible space-like Cauchy hypersurfaces N\ C (MHD) t , N-2 C (MHD) t , 
ti 7^ t 2 , such that the following properties are satisfied: 

(a) N t C (+> +)Yi; 

(b) N 2 C (MffP) ; 

there exists a crossing nuclear critical zone solution V C (MHD), such that dV = 
N\[J P IJ-/V2, where P is a time-like admissible 3-dimensional integral manifold of 
(MHD). In general such a solution is a singular solution. Therefore, under above 

condition of admissibility, one has algebraic classes in fij^ -. Furthermore, V 

is represented by a smooth integral manifold at Y , with respect to the split given in 

(87), if the tangent space T (+) V\ Y = TZ\ Y C E 2 (MHD)\ Y , where E 2 (M HD) is 
the Cartan distribution of (MHD). In fact the dimension of the Cartan distribution 
E2 of (MHD) in the points q G (+.+,+) (MHD) is higher than in the points q G 
( + + )Yi). This can be seen by direct computation. Let us denote by 

,01 _1_ 7.3 p>„,hfa 



0<|j8|<1 

the generic vector field of the Cartan distribution E2(W), where 



{2/ J '}i< J < 1 6 = K,P,^ J B\H\/ < ,P,M 
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are the vertical coordinates of the fiber bundle tt : W — > M. The Cartan distribution 
on the boundary r +t+ \Y\ C (MHD) is given by the vector fields (88) such that the 
following equations are satisfied: ^.F 1 = 0, where F 1 are the functions defining 
equation ( ++ )Yi. These are just the functions defining ( + + , + )(M.ffD), i.e., F^ s ' , 
1 < s < 7, but with the condition h = 0. So we get 

(89) {dimE 2 (Am£>) = dim E 2 ((+,+,+) (MBD) ) = 148 > dimE 2 (( +l+ )fi) = 138.} . 

Note that even if the points o/( +i+ )Yi can be considered singular one, with respect 
to the Cartan distribution of (MHD), the embedding (MHD) C (MHD) allows us 
to prolong a solution from ( + + )Yi to ( + _ +}+ )(MHD), according to Theorem 4. 30 
and surgery representation (87). Then such solution V in general bifurcates along 
Y. By summarizing, we can say that the crossing nuclear critical zone solution 
V C (MHD) is represented by an integral manifold V C (MHD) that is smooth 
in a neighborhood of X C Y, iff for its tangent space TV, the following condition 
is satisfied: T (+) V\ X = TZ\ X C E 2 (MH~b)\ x ■ 

Now, in order to surgery a smooth solution V' C (+, +, +)(MHD), passing through 
a compact smooth space-like 3- dimensional manifold Ni, with X\ = dNi orientable, 
with suitable smooth solutions of r +<+ \Yi it is enough that the following conditions 
should be satisfied: 

( lim^Q V = {+t+) V C ( + , +, +) (MHD) 

(90) <^ 

[ T (+!+) V C E 2 ((+, +, +) (MHD) ) Tl=0 . 

Since ( +-+ )V is a smooth solution we can prolong it to 00, obtaining a solution 
C (+.+ t +)(MHD)^^y There it identifies an horizontal A-plane H C 
E2 ((+_+.+) (MHD)^ ), contained in the Cartan distribution of ( + + + ^(MHD)^ . 
Then H identifies also an horizontal A-plane, that we continue to denote with H, 
in the Cartan distribution of (MHD) +OQ , that on ((+,+) Yi) +rxi , coincides with the 
Cartan distribution of this last equation. Then we can smoothly prolong ^ + ^V^°°^ 
into ((+,+) Yij+oa, identifying there a smooth solution ( + ) V^°°^ . The projection of 
this algebraic singular solution on (+,+) Y\ 1J ^ + + + ^(MHD) identifies a smooth so- 
lution V that has the split surgery property (87). Therefore, it is enough to prove 
that solutions with the property (90) exist in ^ + +t+ ^(MHD). Now, we can see 
that a Cartan vector field £ of ^ + + + -j(MHD) is given in (88) and subject to the 
condition <^.F^ S > = 0. Then conditions (90) are satisfied iff X a h a = 0. Since this 
condition can be satisfied for suitable functions X a on ^ + + + ^(MHD), we recognize 
that in the set of smooth solutions of ( + _ +}+ )(MHD) exist solutions that smoothly 
surgery with smooth solutions of r+ + \Y\. Such solutions are not finite times sta- 
ble in (MHD since there the symbol is not trivial in each components ( +!+ )Yi and 
( + +-+ ^(MHD). However, by using their formal integrability and complete integra- 

bility properties, we can state that in ((MHD) +OQ such solutions are finite times 
stable. So (MHD) is the stable extended crystal singular PDE associated to 

(aThd. 
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APPENDIX A - THE AFFINE CRYSTALLOGRAPHIC GROUP 
TYPES [G(3)] AND [G(2)] 



Table 5. The 32 crystallographic point groups of the space-group 
types [G(3)]. 



Type 


Schoenflies Symbol 


International Symbol 


nonaxial (2) 


C-i — S2, C s — Cih 


I, m 


cyclic (5) 


C n , n = 1,2,3,4,6 


n,n= 1,2,3,4,6 


cyclic with horizontal planes (4) 


C nh , n = 2,3,4,6 


2/m, 6, 4/m, 6/m 


cyclic with vertical planes (4) 


C nv , n = 2,3,4,6 


mm2, 3m, 4mm, 6mm 


dihedral (4) 


D n , ra = 2,3,4,6 


222, 32, 422, 622 


dihedral with horizontal planes (4) 


D nh , n = 2,3,4,6 


mmm, 6m2, 4/mmm, 6/mmm 


dihedral with planes between axes (2) 


D nd , n = 2,3 


42m, 3m 


improper rotation (2) 


Sn, n =4,6 


4, 3 


cubic groups (5) 


T,T h> T d} 0, h 


23, m3, 43m, 432, m3m 



n (n — 1, 2, 3, 4, 6): rotations of 2-k/ti about a symmetry axis. 

n (n — 1,2,3,4,6): rotation n composed with inversion about symmetry centre; m — 2. 

Cn — 7-n — {0, 1, 2, 3, ■ ■ ■ ,n — 1} (cyclic abelian groups); Cj — S2; C s — m; C^i — Sq 

D n - group with an 71-fold axis plus a two-fold axis perpendicular to that axis. 

D n is a non abelian group for n > 2. The group order is 2n. 

D n h: D n with a mirror plane symmetry perpendicular to the n-fold axis. 

D nv - D n with mirror plane symmetries parallel to the n-fold axis. 

O: symmetry group of the octahedron. The group order is 24. One has the isomorphism O — T d . 

0^'. O with improper operations (those that change orientation). The group order is 48. 

T: symmetry group of the tetrahedron, isomorphic to the alternating group A4. The group order is 12. 

The group order is 12. 
Tfi- T with improper operations. Non abelian group of order 24. 
T} r : T with the addition of an inversion. 

A n : group of even permutations on a set of length n. The group order is ~. 



Table 6. Diehdral groups L^m, ro > 1. 



m 


£>2m 


1 


D 2 et Z 2 


m > 2 


D 2m = %m X Z2 

The generator of %i acts on Z m as multiplication by —1 


2 


Di = Z2 x Z2 (Klein four-group) 



Warn that in crystallography diehdral groups are usually denoted by D m — £^2m> m > 2. 



This is just the notation used in Tab. 5. 



Table 7. The 230 affine crystallographic space-group types [G(3)]. 



Syngonies (7) 


Geometric classes (32=Point Groups) 


Bravais types (14) 


Triclinic 


C; = S 2 (1), Ci(l) 


2P 


Monoclinic 


C 2 (3), C lk (4), C 2h (6) 


8P, 5C 


Orthorhombic 


D 2 (9), C 2 „(22), D 2h (2S) 


30P, 15C, 9J, 5F 


Tetragonal 


Gt(6), S 4 (2), C 4h (6), D 4 (10), C 4 „(12), P> 2<i (12), £> 4h (20) 


49P, 191 


Trigonal 


C«(4), S 6 (2), £> 3 (7), C 3 „(6), D 3d (6) 


18P, 7fl 


Hexagonal 


C 3 (6), C 3h (l), C 6fe (2), Dg(6), C 6 „(4), D 3h (4), D 6h (A) 


27P 


Cubic 


T(5), T h (7), 0(8), r d (6), O h (10) 


15P, IIP. 91 



(Crystal class systems —Syngonies): 7. {Triclinic— Anorthic), (Trigonal^Rhombohedral). 

Bravais lattice centering types: P —primitive, / —body, F —face, A/B/C —side, R — rhombohedral. 

(For the body-centerd case (B — I), the infinite translation group is {Z 3 , 7? + {\, ^, §)}■) 
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Table 8. The 17 affine crystallographic space-group types [G(2)]. 



Syngony (5) 


International Symbols 


Geometric classes (10=Point Groups) 


Oblique 


pi, p2 


Zi, Z 2 


Rectangular 


pm, pg, cm, pmm 


Di, D u D l , D 2 


Rhombic 


pmg, pgg, cmm 


D 2 , D 2 , D 2 


Square 


p4, pAm, pAg 




Trigonal 


p3, p3ml, p31m 


Z 3 , D 3 , D 3 


Hexagonal 


p6, pQm 





All planar crystallographic groups G(2) are subgroups of p4m or pQm, or both. {See Appendix D.) 
p —primitive, c —centered, m —mirror plane, 5 —glide reflection. 

A glide reflection is an isometry of the Euclidean plane that combines reflection in a line 
with a translation along that line. 



Table 9. The 4 holohedries (lattice symmetries) in [G(2)]. 



Holohedry 


Lattice system 


■2 


monoclinic=oblique 


2mm 


orthorohmbic— {rectangular, rhombic} 


4mm 


tetr agonal — square 


6mm 


hexagonal— {trigonal, hexagonal} 



Table 10. The 7 holohedries (lattice symmetries) in [G(3)]. 



Holohedry 


Lattice system 


1 


triclinic 


2/m 


monoclinic 


mmm 


orthorohmbic 


4 /mmm 


tetragonal (square) 


3m 


trigonal (rombohedric) 


6 /mmm 


hexagonal 


m3m 


cubic 




FIGURE 4. Octahedron and Tetrahedron 3-chains in R 3 . (The 
tetrahedron is a 3-chain in R 3 , identified with a regular 4- faced 
polyhedron, where each face is an equilateral triangle. It can be 
derived from the octahedron by extending alternate faces until they 
meet.) 
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6. APPENDIX B - THE SUBGROUPS OF THE AFFINE 
CRYSTALLOGRAPHIC POINT-GROUP TYPES [G(3)] 
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Index 






Subgroup 
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Index 
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1 


1 


1 
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Index 
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2 
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1 


-4 


-3 


6 


2 


2/rn 


4 
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2 


2 


2 


3 


3 
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4 
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2 


2 


2/m 


4 


3 


m 


2 


4 
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Subgroup 


Order 


Index 
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6 
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16 
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-42m 
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422 
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Subgroup 


Order 


Index 


23 


12 


1 


3 


3 


4 


222 


4 


3 


2 


2 


6 


1 


1 


12 



Subgroup 


Order 


Index 


432 


24 


1 


23 


12 


2 


32 


6 


4 


422 


8 


3 


4 


4 


6 


3 


3 


8 


222 


4 


6 


2 


2 


12 


1 


1 


24 





Subgroup 


Order 






24 


1 


23 


12 


2 


-3 


6 


4 


3 


3 


8 


mmm 


8 


3 


mm2 


4 


6 


222 


4 


6 


2/rn 


4 


6 


2 


2 


12 


m 


2 


12 


-1 


2 


12 


1 


1 


24 


Oft 


Subgroup 






m-3m 


48 


1 


—43m 


24 


2 


432 


24 


2 


m— 3 


24 


2 


23 


12 


4 


—3m 


12 


4 


3m 


6 


8 


32 


6 


8 


-3 


6 


8 


3 


3 


16 


4/mmm 


16 


3 


-42m 


8 


6 


4mm 


8 


6 


422 


8 


6 


4/m 


8 


6 


-4 


4 


12 


4 


4 


12 


mmm 


8 


6 


222 


4 


12 


2/rn 


4 


12 


2 


2 


24 


m 


2 


24 


-1 


2 


24 


1 


1 


48 





Subgroup 


Order 


Index 


-43m 


24 


1 


23 


12 


2 


3 m 


6 


4 


3 


3 


8 


-42m 


8 


3 


-4 


4 


6 


mm2 


4 


6 


222 


4 


6 


2 


2 


12 


m 


2 


12 


1 


1 


24 



EXTENDED CRYSTAL PDE'S 



o5 



7. APPENDIX C - AMALGAMATED FREE PRODUCTS IN THE 
AFFINE CRYSTALLOGRAPHIC SPACE-GROUP TYPES [G(3)] 



Amalgamat 


ed free 


products in [G(3)] 




5 2 * C Z 2 = 


<x,y> 


: 




: <4,4> 




+2 2 -^2 — 


<2,4,x> 


z e 


★z 3 £>3 = 


<3,6,x> 


D 4 


~kD 2 7)i- 


= <4,4,x> 


_D 2 xZ 


2W D 2 -O4 


= <2,4,I,x> 


D 6 


*D 3 D 6 = 


= <6,6,x> 


D 3 x2 


2~kn s 7D G 


=<2,6,I,x> 



x— reflection over ic-axis; y— reflection over y-axis. 



8. APPENDIX D - THE SUBGROUPS OF THE AFFINE 
CRYSTALLOGRAPHIC SPACE-GROUP TYPES [G(2)] 
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